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ABSTRACT 



An improved formulation for the N(H2)/I(C0) conversion factor or X-factor 
is proposed. The statement that the velocity-integrated radiation temperature 
of the ^^CO J = 1 — >^ line, /(^^CQ), "counts" o ptically thick clumps is quan- 
tified using the formalism of iMartin et al.l (119841 ) for line emission in a clumpy 
cloud. Adopting the simplifying assumptions of thermalized ^^CO J = 1 — * 
line emission and isothermal gas, an effective optical depth, Tg/, is defined as the 
product of the clump filling factor within each velocity interval and the clump 
effective optical depth as a function of the optical depth on the clump's central 
sightline, Tq. The clump effective optical depth is well approximated as a power 
law in Tq with power-law index, e, referred to here as the clump "fluffiness," and 
has values between zero and unity. While the ^^CO J = 1 — line is optically 
thick within each clump (i.e., high tq), it is optically thin "to the clumps" (i.e., 
low Tg/). Thus the dependence of I {CO) on r^f is hnear, resulting in an X-factor 
that depends only on clump properties and not directly on the entire cloud. As- 
suming virialization of the clumps yields an expression for the X-factor whose 
dependence on physical parameters like density and temperature is "softened" 
by power-law indices of less than unity that depend on the fluffiness parameter, 
e. The X-factor provides estimates of gas column density because each sightline 
within the beam has optically thin gas within certain narrow velocity ranges. 
Determining column density from the optically thin gas is straightforward and 
parameters like e then allow extrapolation of the column density of the optically 
thin gas to that of all the gas. Implicit in this formulation is the assumption 
that fluffiness is, on average, constant from one beam to the next. This is also 
required to some extent for density and temperature, but the dependence of the 
X-factor, Xj, on these may be weaker. 
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One important suggestion of this formulation is that virialization of entire 
clouds is iirelevant. The densities required to give reasonable values of are 
consistent with those found in cloud clumps (i.e. ~ 10^ cm~^). Thus virial- 
ization of clumps, rather than of entire clouds, is consistent with the observed 
values of Xf. And even virialization of clumps is not strictly required; only a 
relationship between clump velocity width and column density similar to that 
of virialization can still yield reasonable values of the X-factor. The underlying 
physics is now at the scale of cloud clumps, implying that the X-factor can probe 
sub-cloud structure. 

The proposed formulation makes specific predictions of the dependence of Xj 
on the CO abundance and of the interpretation of line ratios. In particular, the 
^^CO J = 1 ^ 0/^^CO J = 1 ^ line ratio values observed in the Orion clouds 
suggest that e ~ 0.3 ± 0.1. If the majority of the ^^CO J = 1 — > emission 
originates in structures with an density variation, then the constraints on 
e also constrain the ratio of the outer-to- inner radii of the region within 
the clumps. Specifically, this ratio for spherical clumps must be 2 to 9 and 
for cylindrical clumps it must be 4 to 42. This is apparently consistent with 
observations, but higher spatial resolution is necessary to ensure that the observed 
ratios are not just lower limits. This formulation also ties the narrow range of 
the observed values of the ^^CO J = 1 ^ 0/^^CO J = 1 ^ line ratio to the 
relative constancy of the X-factor. 

The properties of real clumps in real molecular clouds can be used to estimate 
the X-factor within these clouds and then be compared with the observationally 
determined X-factor. This yields X-factor values that are within a factor of 2 
of the observed values. This is acceptable for the first attempt, but reducing 
this discrepancy will re quire improving the fo rmulation. While this formulation 
improves upon that of iDickman et al.l (119861 ). it has shortcomings of its own. 
These include uncertainties as to why e seems to be constant from cloud to 
cloud, uncertainties in defining the average clump density and neglecting certain 
complications, such as non-LTE effects, magnetic fields, turbulence, etc. 

Despite these shortcomings, the proposed formulation represents the first ma- 
jor improvement in understanding the X-factor because it is the first formulation 
to include radiative transfer. 



Subject headings: ISM: molecules and dust — Orion 
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1. Introduction 

One of the most basic questions tliat can be asked in any field of researcli is about the 
quantity of the material under investigation. In studies of the interstellar medium (ISM), 
for instance, the amount of gas and dust bears on questions of the physical mechanisms 
that effect and control the ISM and, consequently, how that ISM evolves and can affect the 
evolution of an entire galaxy. In particular, the amount of molecular gas in a cloud, cloud 
complex, spiral arm, or galaxy constrains the number of stars that form and the way that 
they form. The workhorse molecule for e stimating molecular gas masses has been, and still 



is, CO (e.g., see llAU Symp. ^170.1 Il997l . and references therein). Specifically, observations 
of the J = 1 ^ rotational line of the isotopologue, ^^C^^O (just CO for short), permit 
simple, but crude, estimates of the mass of molecular hydrogen in an astronomical source. 
The velocity-integrated radiation temperature, I(CO), often called the integrated intensity, 
is multiplied by a standard conversion factor, N(H2)/I(C0), to yield the molecular hydrogen 
column density, N(H2), which gives the H2 mass of the source after integrating over the 
source's projected area. The most current value of this conversion factor i s about 2 x 



10^° H 2 cm ^ ■ (K ■ km ■ s ^) ^ for the molecular gas in the disk of our Galaxy (jPame et al. 



200l|). 



Why the CO J = 1 — > line should yield an estimate of column density is far from clear. 
Even if this line were optically thin, the conversion factor would depend in a simple way on 
the physical conditions in the molecular gas. Given that CO has many rotational levels and 
that the spacings of these levels (in temperature units) are comparable to the temperatures 
found in molecular gas, the N(H2)/I(C0) conversion factor would depend, at the very least, 
on the gas kinetic temperature. Given that the densities inferred with molecular clouds (of 
at least ~ 10^ H2 cm ^) are comparable to, and not much higher than, the critical densities 
of the observed rotational transitions, then it is obvious that N(H2)/I(C0) should also have 
at least a weak dependence on molecular gas density. Another obvious dependence would be 
on the abundance of CO relative to H2, X(CO). If these physical conditions were known, then 



the m olecular hydrogen column density could be recovered easily (see Appendix A of IWall 



20061 ). Even if these conditions are not known exactly, observations of molecular clouds on 
galactic scales would yield reasonable estimates of molecular gas column densities, because 
reasonable values for the relevant physical parameters are well known and are relatively con- 
stant from source to source. For example, estimating the mass of the molecular medium of an 
entire galaxy using an optically thin molecular line would be, on average, more reliable than 
estimating the mass of a molecular cloud core of unknown physical conditions, because the 
physical conditions in the molecular gas averaged over the scale of a galaxy are less extreme 
and vary much less from galaxy to galaxy than they do from cloud core to cloud core. In 
any event, the conversion from integrated intensity to column density is very straightforward 
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and relatively simple in the opt ically thin case. However, as is well known (e.g., see lEvans 
1980l : lKutnerlll984l : lEvanslll999l ). the CO J = 1 ^ line is optically thick, obfuscating any 
simple explanation as to why it should probe molecular gas column densities. 

Other tracers of molecular gas mass exist, tracers that do not possess the potentially 
serious uncertainties posed by CO J = 1 — » 0. The rotational lines of the isotopologues 
i3q16q g^j^j^ i2q18q Qust ^^QQ g^^^^^ C^^O in short form), for example, can be optically thin 
and, consequently, their integrated intensities have a straightforward relationship with the 
molecular gas column density (provided the physical conditions are known). While poten- 
tially simpler to use for determining column densities, these optically thin lines are nor- 
mally factors of about 3 to 50 weaker than the CO J = 1 — » line (e.g., iKutnerl Il984l : 



Langer fc Penziad Il99d : iNagahama et al.l Il998l : iMaddalena et al.l Il986l ) ; the ^^CO lines are 
better for mapping la rge areas of molecular gas or f or detecting weak sources, s uch as high- 



redshift galaxies (e.g..lBrown fc Vanden Boutlll992l : iBarvainis et al.l 1 199 71 . Il998l : lAUoin et al. 



20001 : ICariUi et al.ll2002l ]all This makes the J = 1 ^ line of CO, and the N(H2)/I(C0) fac- 
tor, more useful or even essential in estimating the total molecular gas mass in some sources, 
resulting in a strong incentive for understanding the N(H2) /I(CO) factor's behavior. 

The usual attempts at accounting for why the N(H2)/1(C0) factor, or X-factor (or 
Xf), is relative l y cons tant on multi-parsec sca les are variation s of the explanation given by 



Dickman et al.l (Il986l ). hereafter DSS86 (e.g., ISakamotd Il996l ) . A summary of the DSS86 
explanation follows. If Tj, is the peak radiation temperature of the CO J = 1 ^ line 
and Av is the appropriately defined velocity width of this line, then I(CO) = Tj^Av. If the 
molecular gas under observation is virialized, then the observed velocity width is related to 
the mass of this gas and, therefore, the gas column density averaged over the solid angle 
subtended by the observed gas. It was then easy to show that N(H2)/I(C0) oc n°-^/Tp. 
The n was the gas density averaged over the virialized volume of gas. DSS86 found that n 
had to be ~ few x 10^ H2 cm~^ to give the observed value of N(H2)/I(C0); therefore it was 
assumed that this volume included entire clouds. Because the CO J = 1 ^ line is optically 
thick, relatively easily thermalized (compared to higher rotational lines of CO), and may 
almost fill the radio telescope's beam at the line peak, it has been assumed that Tp ~ T,,, 



where Tj^ is the gas kinetic temperature, Tj^ (e.g.. iKutner fc Leung|l985l : IWeiss et al.ll200ll ). 

Even if the gas does not fill the beam (a point to which we will return later), we would 

have Tj^ roughly proportional to Tj^ and we would still have N(H2)/I(C0) oc n^'^/T^. The 

basic argument is that the quantity nP'^/T^ does not strongly vary on multi-parsec scales, 

especially due to the weak dependence on density, resulting in a fairly stable value of X. 

Observational evidence does indeed seem to support a roughly constant value of the X- 

f actor to within a fac tor of about 2 for the disk of our Gala xy, where X 
-1 



km ■ s 



[see, e.g. 



2 X 10^" cm"^ ■ (K ■ 

Dame et al.ll200ll : IStrong et al.lll988l . and references therein), although 
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the observations of ISodroski et al.l (Il994l ) and IStrong et al.l (12004) suggest a higher value 



of Xf in the outer disk (a claim that is at odds with ICarpenter et al.lll990l ). The values 
of the X-factor that apply to the disks of other spiral galaxies are often within factors of 



about 3 of that of the Galactic dis 


i X-factor (e.g., 


Young & Scovillel 


1982|; 


Adler et al. 


I992I; 


Guelin et a . 


1995; 


Nakai & Ki 


Lino 


1995; 


Brouillet et al. 


199d; 


Rand et al. 


1999; 


Meier et al. 


2ood 


2001 


; ] 


Boselli et al. 


2002; 


Rosolowski et al. 


boost). 



Nonetheless, there is good observational evidence that the usual value of the X-factor 
does not always apply. In the centers of exte rnal galaxies, the X-factor is factors of about 



5 or n aore lower than the standard valu e (e.g.. iRickard fc Blitzlll985l ; Ilsraellll988l ; IWall et al. 



1993; 


Regan 


2000 


Paglione et al. 


2001 




(ISodroski et al. 


1995; 


Dahmen et al. 


I997I 



20011). as well as in the central region of our own Galaxy 



the galaxy M 3 1 have an X-factor tha t is an order of magnitu de larger than that for the disk 
of our Galaxy (jSofue fc Yoshidalll993l ; iLoinard fc Allenlll998l ). The high-latitude translucent 
clou ds in our Galaxy sho w X-factor variations with a total range of an order of magnitude 



see 



Magnani et al.lll998l ). In some infrared luminous galaxies there is evi dence that the X - 



factor can be roughly an order of magnitude lower than the standard value (lYao et al.ll2003l ). 
In irregular galaxies, the X-factor can be more than an order of magnitude h igher than the 
standard value (e.g.. 



1997; Fu 



mi et al 



found in 



Israeli Il988l ; lOettmar fc Heithausenlll989l ; Ilsraellll997l lal; iMadden et aL 



19991). A deta. iled discu ssion of the shortcomings of the X -factor can be 



Malonev fc Blackl (|l988l )(also see llsraell [l988l : IXutner fc Leunglll985l ). 



Hence, any complete explanation or theory of the X-factor must allow for and account 
for inferred variations of Xf in some cases and, at the same time, relative stability of Xj in 
other cases. This is a difficult balancing act, but achieving such a theory is important for the 
very basic reason that scientific tools must be thoroughly und erstood. While there are a few 
articles that explore the und erlying physics of Xf (e.g., DSS86 lMaloney fc Blacklll988l ; llsrael 



19881 : iKutner fc Leung||l985l ). there are literally hundreds of articles in the literature that use 



or m ention the X-fact or without a detailed examination of its physical properties (including 



even 



Wall et al.l Il996l ). This is in stark contrast to the situatio n with supernoyae type la 



(SNe la), for example. SNe la can be used as standard candles (jPhillips 



1995 



19961; 



Hanauv et al.lll996l Jal) and can constrain cosmological models (iReiss et al. 



1993 



Reiss et al. 



1998 



Perlmutter et al.lll999l ). While these supernovae are often used as standard ca ndles, there 
is also much theoretical and observational work to understand SNe la (e.g. 



2006 



Sauer et al. 



Wandl2007l:lReiss fc Livid I2OO6I : iGarg et all 120071 : 1 James et al.ll2006l : iNeill et al 



2006 



Borkowski et al.ll2006l ). Such work w ill elucidate why S Ne la are standard candles, or why 



they might not be in some cases (see lHowell et al.ll2006l ). In a sense, the X-factor is almost 
a standard candle that relates surface brightness to surface density (i.e., I{CO) to N{H2)). 
And yet there is comparatively little effort to shed light on the physics underpinning the 
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X-f actor. 

Furthermore, there are a number of more specific reasons for understanding Xf-. 



Reliable molecular gas mass estimates. Among other things, this can refine our under- 
standing of star formation yields (i.e., what fraction of the gas goes into star formation). 

Improved knowledge of molecular cloud structure and physics. Better comprehension 
of the X-factor can constrain estimates of molecular cloud physical parameters and can 
supply new insights into star formation processes. 

Improved radiative transfer theory. While it is unlikely that there would be funda- 
mental improvements of radiative transfer theory, such improvements are still possible, 
benefiting astrophysical theory in general. 



These represent possible long-term goals of research into the physics underlying the X-factor. 



The goals of the current paper are considerably more modest: addressing the deficiencies 
of the DSS86 explanation of the X-factor, improving upon this explanation, and examining 
a few consequences of the formulation proposed here. Improvements are necessary because 
DSS86 has the following problems: 

1. No treatment of radiative transfer. This is a fundamental problem with DSS86. At 
first glance, it might seem superfluous to treat radiative transfer in the optically thick 
case. However, if we consider a clumpy medium, where the clumps can have optically 
thin edges and optically thin frequencies in their line profiles, then treating radiative 
transfer is essential for understanding the X-factor. In particular, the optically thin 
limit of CO J = 1 ^ must also he included. Any complete treatment must include 
the optically thin case, whether this case is observed in nature or not. This case cannot 
be included easily in the DSS86 explanation because it includes the virial theorem 
withozit including radiative transfer — virialization by itself says nothing about the 
optical depth of the emission. 



2. Sensitivity to and n(H2). As discussed in IWalll ( l2006l ). I(CO) and the X-factor 
estimate the molecular hydrogen column densities to within factors of about 2 of the 
values determined from optica lly thin trac ers for the majority of positions in the Orion 
clouds. And yet we know from I Wall! (120061 ) (looking at two-component model tempera- 
tures), the range of gas kinetic temperatures can be an order of magnitude. In general, 
we know that molecular cloud kinetic temperatu res and densities have a full range of 



an order of magnitude on raulti-parsec scales ( cf. ISanders et al.lll985l : ISakamoto et al. 



1994 : iHelfer fc Blitzerl 119971 : iPlume et al.l l2000l ) . Since the X-factor supposedly varies 
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as n /Tj^, the temperature and density variations can each change Xf by factors of 
3 to 10 (unless n"'^ were to vary hke Tj^, but this is unhkely to be true in general, 
especially if there is pressure equihbrium) . Thus the Xf of DSS86 is too sensitive to 
the density and kinetic temperature. Having a weaker dependence of Xf on n and 
Tj^, like X oc(n/Tj^)° '^, would resolve this sensitivity problem; variations of an order 
of magnitude in either n or Tj^ would allow X to vary by less than a factor of 2. 

3. Virialization of entire clouds. DSS86 require low densities (i.e. n(H2) ~/ewx 10^ cm~^) 
to obtain the observed value of the X-factor. Given that the critical density of the 
CO J = 1 — > transition is ~ 3 x 10^ H2 cm~^, the densities of the CO-emitting 
structures are abou t an order of magnitude higher (also see the average densities of the 
filaments found by Nagahama et al. llQQsl ). The low density required by DSS86 may 
represent the density averaged over an entire cloud. For the Orion A and B clouds, 
this volume-averaged density is between about 200 and 600iJ2cm~^, depending on 
the precise assumptions used. Therefore, DSS86 are assuming that entire clouds are 
stable and virialized. Some evidence sugge sts that larger molecular clouds are indeed 



virialized (i.e. for masses ^ 10"^ M©, e.g., see lHeyer et al. 



2001: Simon et al. 



other evidence suggest s that many molecular clouds may not be (jPringle et al.l 12001 



Clark fc Bon"^ 120041: 



yazguez-Semadeni et al.l |2007(@I ). especially in extragalactic 



2001), while 



systems (e.g., see Israel 200ol ). Also, the X-factor seems to yield reasonable column 



density estimates for gas on scales smaller than entire clouds (e.g., IWallll2006l ). 



4. Stronger dependence of peak on N(H2) than of Av on N(H2) is not explained. 
DSS86 require that the observed velocity width of the line depends on the gas col- 
umn density. However, there is evidence that it is the peak radiation temperature. 



^R' -f-- 

Figure [1] and 


Wall 


20061: 


Hever et al. 


1996 


Ballesteros-Paredes & Mac Low 


2002 


)• 



19961 : IPichardo et al.ll2000l : lOstriker et al 



2001 



The purpose of the current paper is to propose an improved approach for understanding the 
X-factor that will resolve, or at least mitigate, the problems with DSS86. For example, the 
explanation proposed here includes radiative transfer in a clumpy medium and shows how 
the optically thick CO J = 1 — emission of a cloud can be sensitive to the optical depths 
of the individual clumps. As a result, this explanation will permit, in some circumstances, 
a very weak dependence on Tj^ and n(H2). Also, even though we will also use the virial 
theorem (except in one case), we can apply it to scales smaller than entire clouds. And 
the X-factor in the current proposed explanation will lose its dependence on virialization 
in the optically thin case. In addition, the proposed approach will naturally explain the 
dependence of the peak Tj, on N(H2). This improved approach has shortcomings of its own, 
but nevertheless represents the first major improvement in understanding the X-factor since 
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DS S86, hecause it i s the first formulation to include radiative transfer. (Tlie reader may also 



consult IWallll2006d . for a very brief description of the method.) 



2. A Formulation for the X- Factor 

2.1. Radiative Transfer in a Clumpy Cloud 

The X-factor may yield a reasonable estimate of the molecular gas column density, 
because the integrated intensity of the CO J = 1—^0 lin e is essentially counting optically 



thick clumps in the gas in the beam (e.g., see lEvanslll999l ). Two clumps on the same line of 



sight within the beam will be, on average, separated in velocity by more than the velocity 
widths of the individual clumps, thereby allowing the clumps to contribute their intensities 
to separate velocities within the line profile without absorption of the emission from the 
more distant clump. And clumps at the same velocity within the line profile will be, on 
average, at different locations within the beam, their intensities simply added together at 
that velocity within the profile. This explanation does not, by itself, directly relate the 
masses of individual clumps to the observed integrated intensity, because, again, the clumps 
are optically thick in the CO J = 1 — >^ line. Applying only the DSS86 approach to the 
clumps will not work, because, as discussed in the introduction, DSS86 and the observed 
value of the X-factor together require densities an order of magnitude lower than is found 
in the clumps of real clouds. The DSS86 derivation of the X-factor depends on the beam- 
averaged column density, N, and the observed velocity width, Av. We need a treatment of 
the problem in which the beam-averaged quantities, N and Av, are cancelled out in favor 
of the corresponding quantities for an individual clump, i.e., Nc and Av^. And we need a 
treatment of the radiative transfer in a clumpy medium. 



Martin et al.l (119841 ) (hereafter MSH84) developed a method for describing radiative 
transfer through a clumpy medium in a highly simplified case: they assumed that each 
clump was homogeneous and in LTE. For additional simplicity, they also assumed that 
the clumps were identical, although they pointed out that their method could be easily 
generalized to clumps with a spectrum of properties (see the Appendix of MSH84). The 
assumption of LTE was necessary because the implicit assumption is that the excitation 
temperature of the transition is constant throughout each clump, which is easily attained 
if the density is high enough for LTE. If the density is not high enough for LTE, then 
the populations of the rotational levels are affected by the ambient radiation field at the 
frequency of the line: all other things being equal, the molecules on the surface of a clump 
will be less excited than those in the clump center, because the former only see the radiation 
from roughly 2tt steradians of solid angle, while the latter see it from the full An steradians. 



- 9 - 



This results in a spatial gradient in the line's excitation temperature. If the density is high 
enough for the transition to be in LTE, then this excitation temperature will be equal to 
the kinetic temperature of the gas throughout the clump. Since the clump is assumed to be 
homogeneous, this kinetic temperature is constant throughout the clump, thereby ensuring 
that the line's excitation temperature is also constant. These assumptions are particularly 
appropriate for the CO J = 1 — line: because of its high optical depth (i.e., r ~ few) 
and low critical density (i.e., n^^.^ ~ 3 x 10^ cm~^), this line is largely thermalized (i.e., close 
to LTE). Hence the method of MSH84 is appropriate here. 

MSH84 used a statistical approach to find the appropriately averaged optical depth 
on a sightline through a cloud with clumps in a vacuum. The effective optical depth on 
a given sightline was expressed in terms of the individual clump opacities and the mean 
number of clumps on a sightline with velocities within a clump's velocity width for the 
case of identical clumps. This effective optical depth is the expectation value of the total 
optical depth of the clumps on a given sighthne, considering the probability of a given 
clump impact parameter with respect to the sightline (i.e., perpendicular displacement of 
the clump center from the sightline) and, accordingly, of a given line-of-sight optical depth 
through each clump. Computing this expectation value then depends on an average opacity 
over all impact parameters for each clump, which is the appropriately determined average 
opacity over the clump's surface area projected in the sightline's direction. The appropriate 
average of the optical depth is determined from the average over values of [1 — exp(— r)] 
and not over r itself because the observed emission depends directly on the former and only 
indirectly on the latter — and the relationship between the two is non-linear. Accordingly, 
the effective optical depth, Xe/, is given by 1 — exp{—Tef) = (1 — exp(— r)), which implies 
exp{—Tef) = {exp{—T)), where () indicates expectation value. In the approximation of the 
spectral line width, Av, being much larger than the velocity width of an individual clump, 
^Vc, Tef can be expressed as the product of the number of clumps per clump velocity width on 
a sightline and the effective optical depth of an individual clump. If is the beam-averaged 
gas column density and Nc is the gas column density averaged over the projected area of a 
single clump, then (A^/Nc)(Avc/Av) is the number of clumps per sightline averaged over the 
beam per clump velocity width at the line central velocity. (Note that MSH84 defined as 
the number of clumps per unit projected cloud area rather than the column density of gas. 
Note also that, below, Nc is actually the column density on the central sightline through the 
clump, but this change in definition accords with the definition of the clump average optical 
depth; see Appendix |A] for details.) If tq is the optical depth on a sightline through the 
center of a single clump, then, following MSH84, A{tq) is the clump effective optical depth. 
(Note that MSH84 called A{tq) the effective optically thick area of the clump. Even though 
their term is more accurate, the simpler "clump effective optical depth" is adopted here.) 
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Consequently, 



A(ro) exp 



2Av2 



(1) 



where Vz is the velocity component along the sightline and where a Gaussian line profile has 
been assumed. If T{x,y) is the clump optical depth on the sightline at position (x,y) with 
respect to a sightline through the clump center, then the clump effective optical depth is 
given by 



A(ro) 



1 



^ A 



vc a^,, 



dv / dx / dy < 1 — exp 



-r(x, y) exp 



2 A^ 



(2) 



where a^^.^. is clump's effective projected area defined in terms of its optical depth: 



a,ff = ^ I dx I dy T{x,y) 



(3) 



The To is simply t{x = 0,y = 0), the optical depth through the clump's center and at 
the center of the clump's velocity profile. The integrals are over the projected area of the 
clump and over the clump's velocity profile. For more details, see MSH84 and Appendix lAl 
Figure [2] shows the variation of A(ro) as a function of tq for two types of clumps: cylindrical 
(seen orthogonally to the axis of symmetry) and spherical. 



The observed line radiation temperature, Tj^, is then related to r^f by the usual expres- 



sion 



where 



Tjz/) = X(TJ [l-exp(-ref)] 



hu r 



/ hz/ 

exp 



VkT, 



-1 



/ hz/ 
exp - 

Vi 



,kT. 



- 1 



-1 



(4) 



(5) 



Tj^ and Tg,., are the gas kinetic and cosmic background temperatures, respectively. As stated 
earlier, LTE is assumed for the emission of the spectral line at frequency, u. The J'^{T^) 
is the source function in temperature units and is the correction for the cosmic microwave 
background emission and for the failure of the Ray leigh- Jeans approximation. Of course 
when Tej ^ 1, we have the simplified form of equation (jl]): 



Tjz/)= J,(Tjref 



(6) 



This is often called the "optically thin limit" for the equation of radiative transfer. A very 
interesting and important point here is that the effective optical depth is in the optically 
thin limit even though the individual clumps can still be quite optically thick. And this, of 
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course, will provide a partial explanation for the X-factor. Substituting equation ([T]) into 
equation ([6]) yields 



As mentioned previously, the quantity (A^/Nc)(Avc/Av) is the number of clumps per sight- 
line averaged over the beam within a clump velocity width at the line central velocity. (Note 
that this is not exactly correct. See the last paragraph of Appendix |A] for an explana- 
tion.) This quantity can be much less than unity, thereby permitting Xe/ <^ 1 even for 
A(ro) ^ 1. (When (iV/Nc)(Avc/Av) < 1, it is similar to the geometric area filling factor 
within a narrow velocity interval, although it is not necessarily equivalent.) Since A{tq) 
is roughly equivalent to [1 — exp(— r)] for a single clump, the meaning of expression ([7]) is 
clear: it is the specific intensity of a single clump at velocity Vz — ~ J'^iT^^) [1 — exp(— r)] 
— multiplied by the number of clumps at that velocity within a clump velocity width — 
(A^/Nc)(Avc/Av) exp[— 1^2^/(2 Av^)]. Simply multiplying the intensity of a single clump by 
the number of clumps gives the observed intensity if the clumps are radiatively de-coupled, 
and this is ensured if Tg/ ^ 1. As r^f increases and becomes optically thick, the different 
clumps within each velocity interval start absorbing each other's emission and the radiation 
temperature approaches the source function j7^(Tj^) asymptotically. 

We are now better equipped to understand the behavior of the curves in Figure O A(ro) 
represents the level of emission from a single clump averaged over the clump's projected area. 
When tq <^ 1, A{tq) ~ tq because all lines of sight through the clump and the line's profile 
at all the clump's internal velocities are optically thin. As Tq increases past unity, the line 
profiles on sightlines passing near the clump's center start saturating in their cores and A{tq) 
starts deviating noticeably from the A{tq) = tq line. Nevertheless, A(ro) continues rising with 
increasing tq because sightlines away from the clump's center are still optically thin. Even 
for sightlines near the clump's center, the line remains optically thin at velocities outside the 
line's core. As tq continues increasing, A{tq) deviates further and further from the A{tq) = tq 
line because the area of optically thick emission slowly increases outwards from the sightline 
through the clump center, covering more and more of the clump's projected area. The 
optically thick portion of the line profile on each sightline increases as well. Nevertheless, 
A(ro) continues growing because of those lines of sight and those velocities at which the 
emission is still optically thin. Figure [2] shows two curves: one for a cylindrical clump of gas 
(i.e. a filament) and one for a spherical clump. The cylinder is viewed side-on (i.e. with its 
symmetry axis perpendicular to the sightline) and has length h. If the symmetry axis is the 
then a Gaussian variation of the optical depth with y was adopted: 





8 
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The spherical clump also has Gaussian spatial variation with optical depth, but with radial 
distance, p, from the central sightline through the clump: 

r{x,y) = Toexp ( -TT^'] , (9) 

where p = a/x^ + y"^. This case was also treated by MSH84, and it is included here for 
comparison. (Note that the a^^^ used here corresponds to the of MSH84.) The effective 
optical depth of the spherical clump grows faster with tq for to> 1 than that of the cylindrical 
clump because the former's optically thick area is growing simultaneously in two dimensions, 
whereas the latter 's grows only in one. While A{tq) can grow without bound in these idealized 
cases, Tj^ cannot. Eventually, A(ro) will grow large enough that r^f -C 1 is no longer valid 
and Tj^ asymptotically approaches j7^(Tj^). The growing tq causes this to happen because the 
clumps start crowding each other spatially and in velocity, due to their increasing optically 
thick and their increasingly saturated line profiles. 

The curves of Figure [2] demonstrate that we can represent them as power-laws in tq for 
Tq < 1 or To > 3: 

A(ro)^kAro^ . (10) 

The values of fc^ and e obviously depend on the specific t(x,?/) — the opacity structure of 
the clump, except in the optically thin case. When tq < 1, we have fcyi = 1 and e = 1, 
regardless of the specific variation of T{x,y). A lower value of e, i.e. closer to zero, indicates 
a clump with a better defined outer edge like a hard sphere. Conversely, aside from the 
optically thin case, a higher value of e, i.e. closer to unity, indicates a clump with a more 
tenuous, or fluffier, outer region. Accordingly, e will be called the "fluffiness" of the clump. 



2.2. Relating Clump Velocity Width with Column Density 

DSS86 required virialization in order to relate the line velocity width to the gas column 
density. That is also required here, but it will be combined with the radiative transfer in a 
clumpy cloud discussed in the previous subsection. The virial theorem in its simplest form 
neglects the effects of surface pressure and magnetic fields, yielding 

2T + W = , (11) 

with T as the total internal kinetic energy of the cloud and W as its total internal potential 
energy. Assuming a spherical clump of uniform density gives 
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where Mc and R are the clump mass and radius, respectively. The total kinetic energy is 

T = ^M,Avl, . (13) 

The Av^d is the three-dimensional velocity dispersion of the gas. The exact kind of velocity 
dispersion this is depends on a number of factors, including the radiative transfer through 
the gas (E. Vazquez Semadeni, priv. comm.). Nevertheless, this velocity width most closely 
resembles an rms width. The Av^d is related to the one- dimensional velocity dispersion Avc 
by Atig^ = 3 Av^. If we assume that molecular hydrogen is the only form of hydrogen in the 
gas with number density, ric, then the mass of the spherical clump is 

M, = —n.fim^^R^ . (14) 

mjj2 is the mass of the hydrogen molecule, fi is the helium correction, and R is the clump 
radius. The density ric is related to the column density through the clump center by ric = 
Nc/Lc, where Lc is the path length on the central sightline through the clump is equal to 
2R. Substituting equations f|T3|) . f|T^ . and f|T2|) into f|TT]) gives us 

Ave = KN'/L'^-' (15) 

and 

L. ^ ( 

.15 



kv = ( — G/imjjj) . (16) 



Numerically in cgs units, this is 

k^ = 2.47 X 10"^^ , 

where /i = 1.3 was used. For a more detailed treatment of spherical clumps, see Appendix [Cl 
(Also, see Appendix [D] for a treatment of cylindrical clumps.) 



2.3. Relating Clump Optical Depth with Column Density 

The clump optical depth on the sightline through the clump's center, tq, can be written 
in terms of the column density of CO in level J, Nj: 



in 



exp 



-1 



Tov 27rAvc 



;i7) 



This comes from equation (A9) of lWalll (120061 ) after applying the Boltzmann factor to change 
Nj^i to Nj. The velocity integral was replaced by to^/2ttAvc, where tq is the optical depth 
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at the center of the clump's velocity profile and on the sightline through the clump's center. 
Tj is the energy of the J J—1 transition in units of temperature: i.e., Tj = hVj j_^/k 



with as the frequency of the transition. A 



and A^,^_ 



J are the spontaneous transition 



rate and the wavelength of the transition, respectively. LTE is assumed, so Tj^ applies in 
place of T-j5^(J — >^ J — 1). We can determine the total column density of CO, N{CO), by 
substituting equation (fT7|) into equation (A22) of IWalll (120061 ): 



N{CO) 



(2J+1)A,,_,A3,^_^ 



Q{T^.)exp 



T 



exp\ 



27rAv, 



where Q{Tj^) is the partition function of CO. Setting J to 1 and rearranging for tq results in 



To 



8y2 7r2AveQ(T^: 



exp 



10 



NcX(CO) 



(19) 



The N{CO) was replaced by N^X ( CO), wher e X{CO) is the abundance of CO relative to H2 
The following values are used (see IWallll2006l . and references therein): Aiq 
Tio = 5.54 K, Aio = 0.2601cm, and X{CO) = 8 x 10"^. Accordingly, 



7.19X 10-S-\ 



1.21 X 10- 



-14 



To 



exp 



5.54 



/2^AveQ(T^; 

The above expression can be represented more simply as a power-law in Tj^: 



27rAv, 



NcT"^ 

K 



(20) 



(21) 



The exact values of k^- and 7 depend on the temperature range and can be computed by 
numerically comparing expressions ( 1211) and (l20l) . In the high-temperature limit, however, 
an analytical solution is possible. This limit means that Tj^ '^Tio and Q(Tj^) 2Tj^/rio 
and [1 — exp(— Tio/T^)] Tiq/T^. This results in kr = 1.85 x 10^^^ in cgs units and 7 = 2. 
But we will be interested in the temperature range Tj^ = 10 to 20 K. The necessary numerical 
comparison gives us 



kr = 7.23 X 10~^^ {cgs units) 



and 



7 



1.75 



for that range. This approximation is accurate to within 1-2% on the above specified range. 
Equation ( 12T]) simplifies further by using expression ( ITSj) for Avc and ric = Nc/Lc- 



To 



L,\/2n 



n 



0.5 



K 



(22) 
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This interesting result suggests that the optical depth for this simplified case (i.e., the velocity 
profile of the optical depth is a simple Gaussian) of a virialized clump does not explicitly 
depend on the sightline pathlength nor the vel ocity widt h, but on their ratio. This is 
related to the Sobolev approximation (e.g., see IShul Il99ll ) in which the optical depth is 



^e.g., see 

dependent on the velocity gradient within a given region and not explicitly on the region's 
size. The pathlength-to- velocity-width ratio (Lc/Avc) in a virialized clump is determined by 
the average density, the spatial variation of the density, and the geometry. Therefore, the 
optical depth depends on those things and the gas temperature, but with no dependence on 
the clump size or velocity width (at least for this simplified case). 



2.4. The X-Factor 



Understanding how to combine the results of the previous subsections to derive an ex- 
pression for the X-factor requires examining the observational data t hat inspired the current 
paper in the first place. Figure [3] shows the Orion data discussed in [Wall (120061 ): the peak 
radiation temperature of the ^^CO J = 1 ^ line (i.e., Tj,) for various positions in the 
Orion clouds normalized to the source function at each position (i.e., JTu(T^k)) versus the gas 
column density (i.e., N(H2)) as determined from ^^CO J = 1 — 0. The plots demonstrate 
a clear correlation between Tj^/j7i/(Tj^) and N(H2). The Spearman rank-order correlation 
test indicates that the correlation exists at better than the 99.99% confidence level. (In fact, 
the confidence level of the null proposition of no correlation is zero to within the machine 
precision — 10~^^.) This is more than just the expected correlation between the J = 1 — 
lines ofJfCO and ^'^CO, because the j7v(Tj^) is determined from the dust temperature (see 
Wallll2006l . for details). This suggests that the dust temperature really is a reliable measure 
of the kinet ic ternpera ture of the molecular gas, at least for the Orion clouds on the scales of 



parsecs (see lWallll2006l JalJbl. for more discussion of this). One way of explaining the correlation 
visible in Figure [3] is that the area filling factor of the clump in each clump velocity interval 
is less than unity. A rising beam-averaged column density, N, could mean that the filling 
factor is rising as more and more clumps fill the beam in each velocity interval. Eventually 
the clumps start crowding each other within the beam and within the line velocity profile 
and the T^/ J',^(T^) ratio starts to saturate and asymptotically approaches unity. (Another 
possibility is that rises because of rising Nc within each clump, thereby increasing the 
clumps' optical depths. This would also produce the observed saturation effect without in- 
creasing the number of clumps within each velocity interval.) Obviously, the goal here is to 
be very specific about the relationship between T^/ J'^iT^) and N(H2). There is sufficient 
scatter and uncertainty in the data that it is not possible to rule out a priori a number of 
such relationships. 
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Nevertheless, from simple radiative transfer theory, we know that the specific intensity 
of a source normalized to its source function, usually written J,^/ S,^, will vary like 1 — exp(— r) 
when plotted against the optical depth through the source, r. Given that the column density, 
N, is proportional to r for constant kinetic temperature and density, the data in Figure [3] 
mimic a curve with the form 1 — exp{—aN) (see the plotted curves), where the aN probably 
represents some kind of optical depth. The majority of the data points are on the roughly 
linearly rising portion of the curve. This represents the optically thin region of the curve, 
but the ^^CO J = 1 — s> zs known to be optically thick from comparisons with the optically 
thin isotopologue ^^CO. Therefore, a clue to understanding the X-f actor is realizing that 
CO J = 1 ^ emission behaves like it is optically thin, despite being optically thick. This 
apparent contradiction is resolved when we consider the effective optical depth as described 
previously. While the individual clumps are themselves optically thick in CO J = 1 — > 0, the 
cloud is optically thin "to the clumps." In other words, the emission from every clump in the 
telescope's beam through the cloud reaches the observer. An analogy would be observing 
the HI 21-cm line from an atomic cloud. In this case, the cloud is optically thin "to the 
atoms" in the sense that the emission from every atom in the telescope's beam through the 
cloud reaches the observer. And since every hydrogen atom is nearly identical in its 21-cm 
line emission properties, the conversion from I(HI) to N(HI) is physically straightforward 
and undisputed. For converting from I(CO) to N(H2), assuming absolutely identical clumps 
would give a constant value of the X-factor, relatable to the clumps' properties. However 



assuming identical clumps contradicts observational evidence (see, e.g.. lTachihara et al.ll2000 



Nagahama et al.l Il998l : iKawamura et al.l Il998l : lOnishi et al.l Il996l ) . But we need not restrict 



the clump properties so severely to explain the X-factor. All we need is to have the clumps 
similar on average from one beam to the next for the X-factor to stay relatively constant. 
And, of course, we must also allow the X-factor to vary in some cases (see Introduction); the 
clumps' average properties must vary from the "norm" in some clouds and locations. 

We now need to quantify this picture, so that we might better understand it and its 
limitations. As Figure [3] clearly shows, the T^/J^{TJ oc N for N<1 to 2 x lO"^"^ H2 cm-"^ . 
This is in the Tg/ -C 1 limit, so equation ([7]) applies and it has the desired proportionality. Of 
course, this proportionality is only visible if the clump properties — N^, Av^, and A(ro) — 
and the observed line width, Av, do not vary strongly with A^. In fact, the scatter visible in 
the plots of Figure ([3]) is probably due to variations in all four of these quantities. Integrating 
equation ([7]) over velocity, Vz, gives 

I (CO) = V2^Tj, (0)Av 

= v^J,(Tjref(0)Av . (23) 

Tj^(O) is the radiation temperature of the CO J = 1 ^ line at f ^ = and is also the peak 
radiation temperature of this line. Similarly, Tej(O) is the effective optical depth at = 0. 
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The X-factor is then given by 



X 



f 



-1 

T-l 



27rX(Tjref(0)AvN-i . (24) 



If we now substitute equation ([T]) evaluated at f ^ = into the above, then 

Xf = 



27rX(TjA(ro)Av,N-^J . (25) 

The important thing to notice here is that the directly observed quantities, N and Av, have 
been replaced by the corresponding clump properties, Nc and Avc. In fact, all the parameters 
in expression (125|) are clump parameters, as desired. It is convenient to define 

Ct = 7,1^ , . (26) 

We now use the approximation Jy(!T^ ~ Tj^ — 3.4 K for Tj^> 10 K and the frequency of the 
^^CO J = 1 ^ line, v = 115.271 GHz. This is good to within 0.4% of T^, (and within 0.6% 
of ^/^^(Tj^)). Consequently, 

Ct = — , (27) 

Tj,-3.4K ' ^ ^ 

which approaches unity as Tj^ grows large. Now we substitute the results of the previous 
subsections into equation (125|) : equation (fTOj) for A{tq), (l22l) for tq, (fT5|) for Avc, and Tj^/Cy 
for j75/(Tj^). Except in the case of the end-on cylinder, where we defined the relationship 
between Nc and ric differently, we also use ric = Nc/Lc. These substitutions yield 

Xf = (27r)5(^-i) Ct k'/ '^T^ ^1^'"'^ • (28) 

The expression (l28ll and its variants (e.g., see Appendices) will be examined in detail. It 
should be mentioned that expression fl28|) is more general than for just a clumpy medium 
and can also apply to a uniform-density cloud (see Appendix [B|) . 

The above formulation for Xf obviously accomplishes the goal of insensitivity to the 
parameters T^^ and ric that we have sought for Xf. The fluffiness parameter, e, is in the 
range to 1; any value in that range that is greater than will confer a greater zrisensitivity 
than occurs for the DSS86 explanation. A particularly interesting example is that value of 
e for which 7e — 1 = 0. In the high temperature limit, Ct 1 and 7 — > 2 and, if e = 0.5, 
then Xf has no dependence on temperature. (Notice that the density dependence is also very 
weak in this case: Xf oc n°'^^.) Given that the CO J = 1 — > line is optically thick, this is 
counterintuitive; raising the temperature by some factor should simply increase I {CO) by 
the same factor (in the high-Tj^ limit), thereby decreasing Xf by that factor. That is not the 
case here. Here we are dealing with a clumpy medium where the optical depth varies across 
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the projected area of each clump. There will always be some sightlines through a clump that 
will still be optically thin. There will also be some velocities in the clump's spectral line 
profile where the line emission is still optically thin. As T^^ increases, Tq goes like T~^, so 
that A(ro) goes like T~^ (see equations and [TUl) . But J'l^iT^^) goes like T^ (in this high-Tj^ 
limit), meaning that the observed Tj^ stays constant. The effect of the increasing kinetic 
temperature of the gas is cancelled by the decreasing effective optical depth of the clumps. 
Another way of saying this is that the effect of the rising temperature is cancelled by the 
shrinking effective optically thick areas of the clumps; the filling factors of the clumps decline 
as the temperature rises. (Note that this special case also occurs for lower temperatures. 
For Tj^ = 10 to 20 K, for example, Ct oc T"*^'^^ and 7 = 1.75. The value of e for which Xf is 
independent of temperature would be 0.75.) Therefore, despite the optical thickness of the 
spectral line in the emitting clumps, changing the optical depths of the individual clumps will 
still have an appreciable effect on the line strength. And this will reduce the dependence of 
the X-factor on the temperature and the density of the gas within the clumps. 

In general, the X-factor provides estimates of gas column density because each sight- 
line within the beam has some optically thin gas within certain narrow velocity ranges. 
Parameters like e then allow extrapolation from the optically thin gas to all the gas. 

One problem with the above analysis is that it assumes clumps of homogeneous density. 
This is inconsistent with most of the t{x, y) functions that will be discussed in Section [31 
The full analyses of these cases are given in the appendices. 

3. Examination of the Properties of the X-Factor 

In this section we examine the properties of the X-factor as formulated in the previous 
section. 



3.1. Dependence on Clump Type 



Here we examine how the X-factor relates to properties of the individual clumps, such 
as velocity width, optical depth, and mass. Accordingly, a "standard" clump — or, rather, a 
set of standard clumps — with specific input parameters must be adopted: geometry, dimen- 
sio ns, and density. These i nput parameters are based on the filamentary clumps as identified 
by iNagahama et all (119981 ) with their ^^CO J = 1 map of the Orion A molecular cloud. 
(It should be mentioned here that the units responsible for the CO -line and dust-con t inuum 
emission may be subfragments within the filaments identified by iNagahama et al.l (119981 ). 
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The filaments identified by that paper are nonetheless used as a first test of the X-factor's 
properties.) Table 2 of that paper lists the Orion A clumps and their characteristics, includ- 
ing dimensions and masses. Based on this table, the roughly cylindrical clumps have lengths 
ranging from 1.7 to 21 pc with a mean length of 6.2 pc, and diameters ranging from 0.7 to 
3.5 pc with and a mean diameter of 1.8 pc. Using the mass-t o-volume ratios, t he av erage 



densities range from about 300 to 4000 ■ molecules ■ cm~^. iNagahama et al.l (Il998l ) also 
used ^^CO J = 1 — >^ data to estimate the gas kinetic temperature for each filament. Their 
Table 2 shows values ranging from about 10 to 40 K. 

The adopted kinetic temperature and density values for the standard clumps use the 
above numbers in combination with other considerations. The adopted kinetic temperatures, 
for instance, also consider the numbers found for the Galactic disk at many-parsec scales. 
For example, the maximum temperature expect ed for the molecular gas and its dust o n 



such scales in the Galaxy is about 20 K (e.g., see lSanders et al.lll985l : ISodroski et al.lll994j ). 



Therefore, the standard clumps have adopted kinetic temperature s of either 10 K o r 20 K 



and not temperatures as high as the 30 to ~40 K values in Table 2 of lNagahama et al.l (119981 ). 
The adopted densities for the standard clumps are loosely related to the range of densities 
mentioned in the previous paragraph, but also expand the range to better explore the effect 
(or to demonstrate the lack of it) on the derived X-factor: the adopted densities for the 
standard clumps are 200, 2000, and 20000/^2 ■ cm'^. 

The adopted dimensions for the standard clumps have no effect on the derived X-factor 
or clump optical depth (see equations and [22]) , except in the case of a cylindrical filament 
viewed end-on. And even in that case, the derived X-factor still does not depend on the 
individual values of the diameter or length, but on their ratio (see Appendix [D]) . Both the 
X-factor and clump optical depth depend on the clump temperature and average density, but 
not on clump size. Nevertheless, the clump dimensions still affect the derived clump mass 
and clump velocity width. For the spherical clumps, the adopted diameter is 1.8 pc — the 
same as the mean diameter of the observed Orion filaments. For the cylindrical filaments, 
the adopted length and diameter are 6.2 pc x 1.8 pc — again the same as the mean values 
of the observed Orion filaments. 

Given that three kinds of density variations are considered for the spherical clumps 
and one for the cylindrical filamentary clumps, there are twenty-four standard clumps in 
the set. For each type of spherical clump, each with diameter 1.8 pc, there are six Tj^,nc 
combinations: T^^ = 10 K, and = 200, 2000, and 2000 H2 cm'^ and T^^ = 20 K with those 
same three densities. For the single type of cylindrical clump, with length x diameter = 
6.2 pc X 1.8 pc, there are the same six T^,nc combinations as for each type of spherical 
clump. The different kinds of density variations considered for the spherical clumps are as 
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follows: uniform, Gaussian, and squared Lorentzian. For the cylindrical filamentary clump, 
only the Gaussian radial density variation is considered, viewed from the side (perpendicular 
to the symmetry axis) and from the end (along that axis). All these cases are listed in 
Table [H along with two extra cases discussed in the next two subsections. 

The following subsections and the appendices examine these cases in detail, the results 
of which are summarized in Tables [D - O 

3.1.1. Completely Optically Thick Case: Optically Thick Disks with Flat-Topped Velocity 

Profiles 

For this case we will examine face-on optically thick disks with completely fiat-topped 
velocity profiles. This means that the optical depth of the ^^CO J = 1 — line will be equal 
to To for all lines of sight through the clump and at every point in the velocity profile for 
each sightline. A(ro) would then have the thoroughly familiar form, 

A(ro) = 1 - exp(-ro) . (29) 

In the optically thick case, i.e. ro^2, 

kA = 1. and e = 0. . (30) 

With fiuffiness and kA set to zero and unity, respectively, all dependence on optical depth 
parameters kA, e, kr, and 7 disappears from expression fl28|) . resulting in an expression for 
Xf that is devoid of radiative transfer: 

Xf = {27r)-'^ Crk^'T-'nl . (31) 

Numerically, this is 

Xf{cgs) = 1.62 X 10^^ CTT-^nl , 

or 

Xf{X2o) = 1.62CTT~'nl , (32) 

where the X-factor in (1321) is in units of X20 or 10^° H2 ■ molecules ■ cm~^ ■ (K • km • s^^)^^. 
(Because the expressions [3T] and 132] come from [28], we have actually assumed that the clumps 
are spheres rather than disks; the result is almost the same so long as A(ro) is given by [291) 

This of course is the DSS86 result with the n°'^/Tj^ dependence for Xf. The resultant 
numerical values for Xf are listed in the first row of Table [H Notice that these Xf values 
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are higher than those for the other cases. Also, notice that the only reasonable values for 
Xf occur for a density of 2 x 10^ H2 cm~^. Given that this result only occurs for this highly 
contrived case of optically thick disks with flat-topped velocity profiles, the DSS86 result is 
unlikely. Nevertheless, DSS86 represents a useful limit. 

3.1.2. Completely Optically Thin Case 

This case assumes that the line emission is optically thin on all sightlines through the 
clump and at all velocities within the line profile. Consequently, A(ro) = tq and 

kA = 1- and e = 1. . (33) 

Now with fiuffiness set to unity, dependence of Xf on and ric disappears. So expression fl2Sl) 
is free of the dependence on virialization; this is the exact opposite of the optically thick 
case discussed in the previous subsection. Accordingly, 

Xj = Crk-^Tl-^ . (34) 

For 1^ = 10 to 20 K, 

Xf{X2o) = 1.38 X 10^2 Ct TO- ^5 _ (35) 

These values of Xf are listed in the second row of Table [H Notice that the Xf values are 
lower than those for the other cases. Equation [35] can also be written as, 

N{H2-cm-^) = 1.38xW^ CtT^^'^^ I (CO) {K- km- s-^) . (36) 

Therefore, the relationship between column density and integrated intensity in the optically 
thin case is recovered. (This is easily verified by starting with expression [23] for A(ro) = Tq 
and combining this with expression [2T1) 

3.1.3. Spherically Symmetric Clumps 

Here we examine cases examined by MSH84: the hard sphere, the Gaussian sphere, and 
the squared-Lorentzian sphere. These examples do not necessarily represent the real clumps 
in real molecular clouds, but simply represent an interesting exploration of the parameter 
space that determines the Xf values. In all the cases considered in this paper, the clumps 
are isothermal and in LTE, implying that the clump optical depth as a function of impact 
parameter, or projected radius p, is proportional to that for the column density: Tc{p) oc 
Nc{p). See Appendix [O for a detailed treatment of the spherically symmetric examples. As 
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stated at the beginning of Section 13. ![ the spherical clumps considered here have radii equal 
to 1.8 pc. 

First, the hard sphere example is examined. This is a uniform-density sphere with a 
well-defined edge at radius, R. Given that the density is uniform, the optical depth profile 
over the projected area of the clump is simply the sightline path-length for impact parameter, 
p: 

V 

i?2 



r{p) = To 



0.5 



(37) 



the same as equation (12) of MSH84, with R in place of r^. Using fl?r|) along with ([T]), ([2]), 
and ([3]) yields the curve depicted with plus signs in Figure 5 of MSH84. Measuring the curve 
for log tq > 0.5 results in kA = 1-7 and e = 0.14. Following the derivation in Subsections 12. 2[ 
12. 3[ and 12.41 or in Appendix [C] yields expression (10131) for Xf. The numerical results for 
this case are listed in Tables [1] to [5] inclusive. The hard sphere is the case closest to the 
completely optically thick case and yields results within factors of a few of those of DSS86 
(see Table [T]). The hard sphere has the advantage of giving reasonable Xj values (i.e. 
~ 2X20 for reasonable densities (i.e. ~ 10'^if2cm~^) and temperatures (i.e. 20 K), but has 
the disadvantage of sensitivity to density and temperature (i.e. Xf oc Ct n^''^^T~^''^^) ■ 

The Gaussian sphere has an optical depth profile given by expression (IGlSp . with t{p) 
and To in place of Nc{p) and uq y/n, respectively. Using this t{p) and the equations (111), (E]), 
and ([3]) yields A(ro) versus tq: the solid curve in Figure 5 of MSH84 and the dashed curve 
in Figure [2] of the current paper. Measurement of this curve gives = 1.6 and e = 0.36 for 
Tq > 3. The derivations in Subsection 10.21 of Appendix O give the numerical results listed in 
Tables [H to [5l As seen in Table [U, the expected Xf values for li = 2 x 10^ H2 cm~^ are within 
about a factor of 2 of the observed value for the Galactic disk. Also, expression (1019P shows 
that Xf oc Ct ii-^'^'^T^^'^^ , which is a greater irjsensitivity to physical conditions than occurs 
for the hard sphere case. 

The squared-Lorentzian sphere's optical depth profile is given by equation fl02ip . with 
a similar expression (i.e. I022p for the column density profile. Following the usual procedure 
gives the A{tq) versus tq curve: the "x" symbols in Figure 5 of MSH84, for which we find 
that kA = 1.5 and e = 0.57. This particular optical depth profile has the advantage of low 
sensitivity to density and temperature (i.e. Xf oc n°'^^T~°-°°^), but the disadvantage of 

/ 



Xf values of roughly factors of 4 too low for a density of ~ 10'^ cm ^ 



In addition to the Xf values (i.e. Table [T]), the tables also have the Avc, tq, A(ro), 
and Mc values (i.e. Tables [2] through [5] inclusive). The Avc and Mc in Tables [2] and O 
respectively, can be compared with the velocity widths and masses of the clumps in the 



Orion A cloud, as listed in Table 2 of iNagahama et al.l (119981 ). Given that spherical clumps 
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with a diameter of 1.8 pc were adopted here, we need the roughly analogous clumps in 



Nagahama et al.l (119981 ) . These are the clumps numbered 14 and 23 in their Table 2, with 



dimensions 2.3 pc x 1.8 pc and 1.7 pc x 1.5 pc, respectively. The observed masses of these 
clumps, 330 Mq and 220 Mq, suggest that they most closely resemble either the hard sphere 
of density 2000 cm~^ or the Gaussian sphere of average density 200 cm~^. However, their 
observed velocity width of 1.1km ■ s~^ suggests that they more closely resemble the hard 
sphere, but with about 2/3 the density, or ^ 1300 cm~^. Here we have also considered that 
the masses given in Table 2 of iNagahama et al.l (119981 ) do not include the correction for 
helium. This lower density would imply Avc ^ 1.2km ■ s^^, close to that observed. This 
implies that these clumps are virialized, as is required for a roughly con stant X-factor to 
apply . However, hard spheres with Tj^ 
1998h and = 1300 cm'^ 



20 K (see column 5 of Table 2 of iNagahama et al. 
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not entirely representative of the rest of the clumps of the Orion clouds. 

Tables [3] and m give the optical depths on the central sightline, tq, and the clump effective 
optical depths, A(ro). The tq values hsted in Table [3] seem to be too high. The observed 
I(^'^CO)/I(^^CO) values for the Orion clouds (see Figure^]) are mostly in the range 0.1 to 0.4, 
implying r(^^CO) of 6 to 24. However, Table [3] lists values that are an order of magnitude 
larger. Nevertheless, given that a roughly constant X-factor only occurs on large scales 
(i.e. many parsecs) and th at the data for which the X-factor was estimated had a spatial 



resolution of about 1° (see IWalll l2006l ). the relevant optical depth to use would be A{tq) 



and not tq. The former is listed in Table H] and has values in the desired range and lower. 
(Although interpretation of the I(^^CO)/I(^^CO) ratio is a little more complicated than in 
the case of homogeneous gas. See Section [3731 for more details.) 



3.1.4- Cylindrically Symmetric Clumps 

Filaments or cylindrical clumps are discussed in detail in Appendix [Dl The behavior 
of these clumps — e.g., the derived tq and Xf values — depend on whether these filaments 
are observed side-on (i.e., perpendicularly to their symmetry axes) or end-on (i.e., parallel 
to these axes). The filament considered here has a Gaussian density variation with distance 
from the central axis, viewed side-on and viewed end-on. As mentioned near the beginning 
of Section Em the adopted dimensions of the filament are 6.2 pc x 1.8 pc. 

For the case of the Gaussian cylinder viewed side-on, the optical depth, r, as a function 
of projected distance from the central axis, y, is similar to equation (ID8|1 . but with T{y) in 
place of Nc{y) and tq in place of Nc{0). Using r(y) with ([T]), ([2]), and yields the solid 
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curve in Figure [2l Measuring the curve for log To > 0.5 gives us fc^ = 1-5 and e = 0.25. Side- 
on cylinders are expected to be less fluffy than spheres with the same optical depth profile, 
given that the former have optical depth variations in only one projected dimension and the 
latter have those variations in two projected dimensions. Here we see that the Gaussian 
side-on cylinder has a fiuffiness that is only about 2/3 that of the Gaussian sphere (i.e. 0.25 
versus 0.36). The derivation in Appendix [Dl Subsection ID. l.lt gives the numerical results in 
Tables [1] to [5l As in the case of the Gaussian sphere, the resultant Xf values for an average 
density of 2 x 10^H2 ■ cm~^ are within a factor of 2 of the value inferred for our Galaxy. Like 
the hard sphere case, this Xf is a little too sensitive to Tj^ and n: Xf oc Ct 



^-0.56 ^0.37 



The Gaussian cylinder viewed end-on is similar to the Gaussian sphere in its optical 
depth profile. However, the central sightline optical depth, tq, depends on the cylinder length- 
to-diameter ratio or aspect ratio, h/di/2- Using the same kA and e as the Gaussian sphere, 
the derivation in Appendix [Dl Subsection ID.2.1t gives the numerical results in Tables [1] 
to [5l The Xf value for an average density of 2 x 10^H2 ■ cm~^ and Tj^ = 20 K is very 
close to that observed for our Galaxy. Also like the Gaussian sphere, this Xf also has the 
advantage of relative insensitivity to Tj^ and n: Xf oc Ct ii°'^^T~°'^^. Unfortunately, these 
advantages are merely lucky coincidences, because the Xf is dependent on the filament's 
aspect ratio {Xf oc which ca n vary by an ord e r of m agnitude or more. For 

example, the filaments of Orion A listed by iNagahama et al.l (119981 ) have aspect ratios that 
vary from about 1 to 10. Another coincidence is that the end-on case implicitly assumes 
that the cylindrical filament is absolutely straight and pointing along the sightline — a low 
probability event. 

Comparison of the results in Tables [1] to [5] with the clump properties listed in Table 2 of 



Nagahama et al.l (119981 ) suggests that the Orion A filaments behave more or less like Gaussian 
cylinders, except for the inferred Xf values. Clumps 12, 22 and 33 have dimensions similar 
to the standard dimensions adopted here: 6.2 pc x 1.8 pc. The average densities of these 
filamentary clumps, n, can be inferred from the column densities divided by the product of 
the observed thicknesses and the for a Gaussian cylinder (see Appendix[D]). These average 
densities and the dimensions can be combine d to give the mass e s and velocity widths. These 



are found to agree with those in Table 2 of iNagahama et al.l (119981 ) to within 7%, except 
for the velocity width of clump 33. The observed value for this clump is about double 
the virialized velocity width, strong evidence th at this clump is not vi rialized. Another 



problem is the inferred Xf values. Table 2 of INagahama et al.l (119981 ) provides enough 
information to roughly estimate the Xf that would correspond to each of the filaments 
listed. Assuming side-on Gaussian cylinders for all 39 filaments listed gives Xf = 0.6 to 
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observations is to consider uniform- density filaments. As seen in Table [T] for the spheres for 
n = 2 X 10^ cm~^, the uniform-density case (i.e., hard sphere) has Xf about a factor of 2 larger 
than for the Gaussian sphere. Accordingly, uniform-density filaments would correspond to 
higher X-factors, closer to that observed for Orion. One slight problem with uniform- density 
filaments is that the inferred virialized velocity widths are lower than those for the Gaussian 
filaments, resulting in a slightly greater number of non-virialized filaments. For uniform- 
density filaments, roughly half have velocity widths within 40% of the virial velocity widths. 
For Gaussian filaments, roughly half have velocity widths within 30% of the virial velocity 
widths. 



3.1.5. Gravitationally Collapsing, Magnetized Filaments 

The examples of clumps that have been examined up till now have all been virialized 
clumps. Strictly speaking, virialization is not required for determining Xj; all that is required 



i s a re lationship between the clump velocity width and the clump's density. iTilley fc Pudritz 



( 120031 ). for example, examine gravitationally collapsing, magnetized filaments. They inves- 
tigate filaments with constant toroidal flux-to-mass ratio and those with constant thermal 
gas pressure to magnetic pressure. In the current paper only the former is considered (see 
Appendix [Dl Subsection ID. 30 . Their equation (35) is rearranged to give the clump velocity 



width in terms of the density; see equation ( lD2ip . iTilley fc Pudritzl (120031 ) find that the 
density goes like r~" at large r, where a = 2 in the case of strong magnetic fields and a = 4 
for weak magnetic fields. Both of these cases are examined here. 

For the case of p(r) oc r^^ (see Appendix |D], Subsubsection ID.3.ip . the total clump 
mass, Mc(ri), depends on the ratio, ri/ro, where ri is the clump's outer radius and Tq is 
the radius of the p = constant core. In contrast, all the previous clump examples had finite 
Mc(cxd). This dependence on ri/ro also extends to A(ro). Numerically integrating ([2]) and 
([3]) shows that e is an increasing function of this ratio, asymptoting out at about 0.9 for 
'"i/'"o~5 X 10^. It is clear then that the fiuffiness of the filament depends not only on the 
density exponent, a, but also on the ratio ri/ro. As was found previously, increased fiuffiness 
had the advantage of reduced sensitivity to Tj^ and li, but the disadvantage of unrealistically 
low Xf. For example, for ti/tq = 1000, we have kA = 1.46, e = 0.69, and Xf oc Ct T°'^^ 
To determine a value for Xf, we must specify a reasonable number for the fragmentation 
wavelength, Xfrag- The simplest way to specify this is to adopt the average cylind er length. 



6.2 pc, that has been used up till now. Equation (35) of ITilley &: Pudrita (120031 ) allows us 



to estimate this number. The observed Avc(FWHM) in most of the filaments in Orion A 
is 1-2 km ■ s^^ or velocity dispersions of 0.4-0.8 km ■ s~^ (after dividing by V8ln2). The 
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conversions from the central density to the average density are found in Appendix [Dl These 
together yield Xjrag = 4.4 to 8.8 pc for a = 2 (and 3.5 to 7.0 pc for a = 4). These numbers 
suggest that choosing Xfrag = 6.2 pc is acceptable. For reasonable temperature and density, 
Tj^ = 20K and n = 2000 ■ cm~'^, Xf = 0.86X20, or a factors of ~5 smaller than that 
for the Galaxy. On the other hand, if ri/ro = 10, we have kA = 1.74, e = 0.27, and 
Xf oc CTT-°-53n0.37_ Tj, = 20K and n = 2000/^2 ■ cm'^, Xf = 1.6 X20, which is 
close to the Galactic value. One problem with the ri/ro = 10 case is that the power- law 
approximation for A{tq) is noticeably poorer than for the othercases studied in the current 
paper. For these other cases, A(ro) ~ kATo is good to within about 13%. For ri/ro = 10, it 
is only good to within about 26%. 

For the case of p(r) oc (see Appendix [D| Subsubsection [D73721) . Mc{oo) is finite and 
the power-law approximation for A{tq) is as good as that for the majority of cases discussed 
in the current paper. Numerical integration of ([2]) and ([3]) yields kA = 1.15 and e = 0.37, 
which results in Xf oc CtT^^'^^ "n.^'^"^- For reasonable temperature and density, Tj^ = 20K 
and li = 2000 H2 ■ cm~^, Xf = 1.9 X20, which is approximately the value for molecular clouds 
in the Galactic disk. 

One complication with the above analysis is that magnetized filaments have velocity 
dispersions that vary spatially wi t hin th e filaments. The central velocity dispersion used in 



equation (35) of iTilley fc Pudrita (120031 ) may di ffer from the obse r ved y elocity dispersion by 



a factor of about 2, as suggested by Figure 4 of iFiege fc Pudrita (120001 ). If so, then Xfrag is 
a factor of two lower than the numbers used above. Of course, if this is true, then the value 
used for k^ must be increased by this factor of 2. Since Xf is a function of the product, 
kv Xfrag, the factor of two decrease in Xfrag is cancelled by this factor of increase in k^, 
leaving Xf unchanged. 

The density dependences found for the collapsing, m agnetized filament are t hose ex- 



pected for an isothermal filament in hydrostatic equilibrium (ITilley fc Pudrital2003l ). lOstriker 



(119641 ) solved for the density distribution in an isothermal filament and found the p{r) oc 
at large r. This solution would also permit p(r) ~ r^^ over a range of intermediate r. 
Since this is similar to the case of the collapsing, magnetized filament, how do the Xf 
values compare? The kpj values were computed in Appendix [Dl Subsection ID.3[ allowing 
determination of k^ for a virialized filament. Using the virialized fc^ values and removing the 
[di/2/ (Xfrag kmax)]^^~^'' factor wiU give the corresponding virialized Xf values. For p(r) oc r~^, 
Xf for the virialized, non-magnetic filament is up to nearly 20% smaller than that for the 
collapsing, magnetized filament. For p(r) oc r~^, Xf for the virialized, non-magnetic filament 
is almost 50% larger than that for the collapsing, magnetized filament. The X-factor values 
derived in the isothermal, virialized case are then easily within a factor of 2 of those for the 
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gravitationally collapsing, magnetized filament. 



3.2. Other Effects on the X-Factor 

3.2.1. Fluffiness 

One advantage of the current, proposed formulation of the X-factor is that Xf can 
be insensitive to the kinetic temperature and density. But this advantage is coupled to an 
important disadvantage: Xf is sensitive to the fluffiness, e. Another way of saying this is that 
the problem of sensitivity to temperature and density has been recast as one of sensitivity 
to another physical parameter. If real clumps in real clouds can indeed be characterized, 
at least approximately, with a parameter e and if the proposed formulation of the X-factor 
is roughly correct, then a relatively constant Xf requires constant e. Understanding the 
X-factor would then provide insights into molecular cloud structure. More specifically, the 
X-factor would itself be a probe of the opacity structure in CO-emitting clumps. 

One obvious way to do this is to plot Xf as a function of e and see what range of e values 
gives realistic Xf values. This must be done for specific and realistic values of Tj^ and li, but 
the parameters fc^, ky, and have implicit dependences on e. These dependences are not 
simple, given that they are affected by many things, such as geometry and density variations. 
However, these dependences are not strong and a crude variation of Xf versus e can be 
plotted if we adopt "typical" values for these quantities: kA = 1.4, k^ = 3.4 x 10"^^ cgs, and 
k^ = 1.5. These values along with li = 2000 H2 cm~^ and T^ = 10 and 20 K were substituted 
into equation (IClip to yield the curves of Figure HI Note that the curve is about 30% lower 
than the level expected for e = and e = 1 (cf. Table [1]). This occurs because k^, k^, 
and especially k^ were fixed. Another point to consider is that the positions of these curves 
also depend on exactly how the "average" densities are defined. Nevertheless, the curves of 
Figure S] still permit approximate estimates of e in molecular clouds. 

Figure H] can estimate the range of e that results in Xf within a factor of 2 of that for 
Galactic disk molecular clouds. Given that Galactic disk molecular clouds have Xf ^ 2 X20 



(iDame et al.ll200ll ). what range of e permits Xf to be in the range 1 to 4 X20? For Tj^ = 20 K 
and li = 2000 cm~'^, e = 0.0 to 0.35. (Taking into account the underestimate of Xf for 
e = 0, this lower limit is more like 0.05.) For Tj^ = lOK and li = 2000 -ff2 cm~'^, e = 0.17 
to 0.46. For Xf = 2X20, e = 0.15 when = 20K and n = 2000i/2cm-3 and e = 0.31 
when Tj^ = lOK and li = 2000if2cm~^. So for the adopted density, the clumps in Galactic 
disk molecular clouds behave most like hard spheres if Tj^ = 20 K and roughly like Gaussian 
spheres or filaments if Tj^ = 10 K. If instead the clumps are filaments with p oc r~^, then the 
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estimated e values constrain the ratio of the maximum-to-minimum radii of the region, 
ti/tq. The e = 0.15 for Tj^ = 20K and li = 2000i72cm~^ requires ti/tq = 2. The e = 0.31 
for Tj^ = lOK and li = 2000i^2cm~^ requires vi/rQ = 16. 

In short, the hkely range of e is subject to a number of assumptions, but is probably 
e = 0.05 to 0.46. One point that should be emphasized is that the curves of Figure HI for 
the given densities and temperatures, are not universal. These curves are dependent on the 
definition of average density (see Section 14.7.11) . 



3.2.2. CO Abundance 



The effect of the CO abundance on the X-factor has been discussed in the literature 



[e.g., see 



Kutner fc Leungj|l985l : iMaloney fc Blacklll988l ). One might expect that Xj would 
depend only weakly on the CO abundance given that th e ^^CO J = 1 — line is optically 



thick. However, as discussed in iMaloney fc Blackl (Il988l ). a reduced CO abundance means 
less self-shielding from the interstellar radiation field and a smaller CO-emitting volume, 
resulting in a lower area filling factor. The current proposed formulation would suggest a 
very specific dependence on the CO abundance, X{CO). From equations f|T9|) and fl2Tl) in 
Section [231 we see that k^. oc X{CO). Equation fl25]) tells us that Xf oc Accordingly, 
Xf oc X{CO)^'^. Given that the most likely values for e are between 0.05 and 0.46 (see 
Section [3.2.11) . this is a weak dependence on the CO abundance. 

However, a more relevant approach is to determine the dependence on the volume- 
averaged abundance. In irregular galaxies it has been found that CO is virtually absent over 
large volumes of molecular clouds, while havin g roughly Ga l actic abundance of CO within 
small central regions within those clouds (see [Israeli [l997al . [2000[ . and references therein). 
Accordingly, an effective CO abundance, X{CO), is defined which is related to the Galactic 
CO abundance, Xg{CO), within a spherical central region of radius, r^o, within a spherical 
clump of radius, ri, by 



X{CO) 



XciCO) { ^ 



(38) 



It is easy to show that the optical depth on the central sightline through the clump, tq, is 
proportional to Tco/ti. From equations ([7]) and ffTOj) we know that. 



Ta(0) OC ^ r, 



(39) 



The quantity NfNc is the number of clumps per sightline averaged over the beam. It had 
been assumed up till this point that the entire volume of each clump was emitting in the 
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CO line. If we now only consider the case where the projected emitting area is less than the 
entire projected area of each clump, then it is easy to show that A^/Nc oc {jcohiY- Putting 
all of this together yields, 

Xf oc X((:70)-3(^+2) . (40) 

Given that e = 0.05 to 0.46, 

Xf oc X(CO)-°-^ *° . (41) 

This is a considerably stronger dependence on the CO abundance than was derived in the 
previous paragraph. In fact, the derived dependence is nearly as strong as that expected 
for optically thin ^^CO J = 1 ^ 0: X/ oc X{COy^. Notice that (gO]), in fact, has this 
dependence in the optically thin case: i.e., when e = 1. 

Filamentary clumps may imply a different dependence of Xf on X{CO) from that 
of spherical clumps. If the filaments are much longer than their diameters (i.e., high as- 
pect ratio), then the dependence of X{CO) on Tco/ti may be two-dimensional rather than 
three-dimensional. This is only an approximation becau se the strong r adiat ion fields or low 



metallicities that can lead to low effective X{CO) (see [Israeli Il997al . l200d . and references 
therein) would reduce the length of a filament's CO-emitting volume as well as its diameter. 
However, for a filament with a high aspect ratio, the fractional change would be larger for the 
diameter than for the length. As long as the effective X{CO) is not extremely low, we can 
roughly assume that X{CO) oc {rco/riY. If the filaments are viewed side-on, iV/Nc oc Tcohi- 
Using a similar approach to that for the spherical clumps. 



Given that e = 0.05 to 0.46, 



Xf oc X(C0)^5(^+i) . (42) 



Xf oc X(CO)~'^-^*°""-' . (43) 



For end-on filaments it is simple to show that, 

Xf oc X{CO)-^ , (44) 

regardless of the value of e. Since filaments, on average, are viewed neither completely side- 
on nor completely end-on, the expected power-law index would be between —0.5 and —1.0. 
Again, notice that gives the dependence expected in the optically thin case for e = 1. 
For (jUl) this optically thin dependence holds regardless of the value of e, given that the 
effective X{CO) is affected only by the projected area of the CO-emitting volume. 
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3.2.3. Effective Optical Depth 

The proposed formulation for the X-factor assumes that the ^^CO J = 1 — > hue is 
optically thin "to the clumps": i.e., Te/ <C 1. But this is not necessarily the case. The 
panels of Figure [3] suggest that saturation starts for N{H2)^ 1 to 2 x 10^^ cm~^. Finding the 
^^CO J = 1 — s> brightness that corresponds to this surface density is not straightforward; 
the inferred Te/ also depends on the source function, J'^(T^). Another consideration is that 
the conditions for saturation that are inferred here really only apply to the Orion clouds 
and not necessarily to other clouds. Still another point is that the peak Tj^ is more relevant 
than the I(CO), because the latter includes the velocity width of the line. The I(CO) of 
an external galaxy, for instance, can have a very large velocity width often dominated by 
large-scale systematic motions, such as the large-scale rotation of the galaxy, rather than the 
smaller scale virialized velocity dispersions within individual clumps or clouds. 

With these caveats in mind, it is at least possible to specify a rough minimum peak Tj^ 
as a necessary (and certainly not sufficient) condition for the start of saturation. Within 
the Orion clouds, this is Tj^^6 to 8K, although there are still positions for which Tj^ ~ 10 
to 12 K and Te/ <^ 1 still applies. So some measure of the source function is necessary to be 
certain that saturation is occurring. 

If saturation is indeed a problem for some sources (or some positions within a source), 
then the actual X-factor necessary for determining the gas surface density could be factors 
of roughly two or more higher than the "standard" value applicable to most other sources 
(or to most other positions within a source). This saturation behavior of the X-factor can 
be understood better by examining either panel of Figure [31 The X-factor of a given point 
on the curve is proportional to the reciprocal of the slope of a line-segment joining the origin 
to that point. As the column density increases, the point moves to the right and the line- 
segment to the origin has an ever decreasing slope; the X-factor increases without bound. 
Therefore, in this optically thick limit, the ^^CO J = 1 ^ line loses its sensitivity to column 
density, as expected. 



3.2.4. Interclump Gas 

Including the effects of a more or less continuous low-density medium between the 
clumps — an interclump gas — in the current treatment of the X-factor is beyond the scope 
of the current paper. Nonetheless, the effects of such gas on Xf can be crudely estimated. 
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From expressions fl23l) . ffTOj) . and ([T]), we have 

/(CO) = v^X(Tj^Av,kAro^ , (45) 

which assumes Xg/ <^ 1, as used throughout most of this paper. Of course, the beam-averaged 
column density, N, is a measure of the mass of CO-emitting gas within the beam. If CO- 
containing gas is added to the beam, then I{CO) increases. How I{CO) increases depends 
on how this gas is added: 

1. If the gas is added in the form of more clumps to each velocity interval, without chang- 
ing the clumps' properties, then I{CO) oc and Xf = N/I{CO) remains constant. 

2. If the added gas increases Nc of each clump, without changing the clumps' dimensions, 
then iV/Nc remains constant and Avc oc N°'^. Given that tq oc Nc/Avc, I{C0) oc 
jY0.5(i+e)^ requiring that Xf oc iV0.5(i-e)_ 

3. If a very distended envelope is added to each clump, then Nc changes very little and 
Avc oc N"'''^. This is a dependence of Avc on because the total mass in the beam 
has increased, which will change the velocity widths of the clumps, even though N^. has 
changed little. This is an inverse dependence of Avc on because the mass is added 
at large distances from each clump center (see more discussion of this below). Again, 
it results that I (CO) oc Aro.5(i+.) Xf oc N^-^^^-'l 

4. There are other cases where Nc is increased only a small amount, but Avc is increased 
by a bit more, resulting in I{CO) oc A^^, where (3 > 1. (Note that if tq < 1, then 
[3 = 1, because e = 1 and tq (Avc/Nc) = 1.) Hence, Xf oc A^-*^"^ and Xf decreases with 
increasing A^. 

One thing to notice about all these cases is that, if the gas is optically thin in the ^^CO J = 
1^0 line (i.e. tq < 1 and e = 1), then /(CO) oc A^ and Xf is constant. Except for point #4 
above, the addition of extra CO-emitting gas will increase, or leave unchanged, the X-factor. 
And point #4 does not necessarily represent a very common case. 

A quantitative estimate of how strongly the interclump gas could affect Xf can come 
from assuming that such gas behaves like an envelope for each clump. The self-potential 
energy of such a clump-envelope system would be 

W = -kwyw , (46) 

Pl/2 

where kw and Mc are the quantities that apply to the clump without the envelope and yw 
is the correction due to the envelope. The internal kinetic energy is 

T=^M,Avla + rM) , (47) 
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where 

ru ^ jI (48) 

is the mass ratio of the envelope to the clump without the envelope. The tm is also 

= Tpvl , (49) 



with 



r, ^ ^ (50) 

PcO 



as the ratio of the central densities of the envelope to bare clump. The is the ratio of 
sizes of envelope to bare clump. Employing the Virial theorem in its simplest form and 
Avj^ = 3 AyI yields 

Avc = kvy°-^ndi/2 , (51) 

with 

Vv = 7— • (52) 

1 + rM 

The central column density of the clump-envelope system, Nc, is 

Nc = 2fcAr?/ArIipi/2 . (53) 

The k^- is the quantity for the bare clump and 

yN = I + TpT^ (54) 
= 1 + rur-^ (55) 

From expressions fl5T]) and fl53l) we see that the and kjq in expression fICllI) for Xf must 
be replaced by k^y^'^ and kjsry^, respectively. The formula for Xf ends up with an extra 
factor of y^f^ in which 

yx = yNy;''-' • (56) 

We are now ready to estimate the effects of the interclump gas on the X-factor. To 
effectively be the interclump gas the envelope must extend far from the clump; i.e., r„ oo. 
Obviously, this yields yw — 1 and y^ ^ {1 + tm)"^- Hence the addition of the envelope 
to the bare clump actually reduces the clump velocity width (see point #3 two paragraphs 
back). Because this limit implies that yN — ^ 1, we have yx — > (1 + tm)^'^- It is clear then 
that the interclump gas has only a small effect on Xf. If the interclump gas mass does not 
dominate the mass of the cloud, then tm < 1 and Xf would increase by less than about 40%. 
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For Tfj not large, these numbers would be smaller. Also, the envelope would increase the 
effective e and /c^ of the clumps, thereby further decreasing the effect of the intcrclump gas 
on the X-factor. (A numerical test of this found that e can increase by about 50%, requiring 
r„ = 10 and the extreme tm = 10. A less extreme tm = 1 produced the usual e but a 25% 
increase in /ca-) 

Knowing the size of the effect on Xf for finite r^- requires knowing the specific density 
variation within the clump and envelope. It is easy to show that for r„ = \ that |/x = (1 + 
Tm^^'^ as it does for r^- oo. (This comes from y^y = (1 + '"m)^ when r„ = 1.) Therefore, for 
To- between 1 and cxo, yx must achieve a minimum or a maximum. (And, since the envelope 
must surround the clump and not vice- versa, must be greater than unity.) Examining 
the cases of a uniform-density spherical clump and envelope and a Gaussian spherical clump 
and envelope shows that both these cases have a minimum for yx in the desired range for 
To-. In the former case, this minimum is greater than unity but less than (1 + rM)°'^- In the 
latter case, the minimum of yx is less than unity but greater than (1 + tm)""'^- Therefore, 
the largest deviations from unity for yx occur for large. Accordingly, only for large is 
the effect on the X-factor maximized; and even in this large-r^ limit the effect is no more 

than 40% if tm < 1- (Note that r^j large means that ^ when rjv/ < 1 and r„ ^ tm 

1 

when tm > i for the uniform-density case and ^ when tm < 1 and ^ r]^ when 
tm > ^ for the Gaussian case.) 

Another point to consider is if the interclump gas has a different temperature from that 
of the clumps. Since this interclump gas is at a lower density than that of the clumps and, if 
there is even very crude pressure equilibrium between the clumps and intcrclump medium, 
then Tj^ (intcrclump) > Tj^ (clump). The warmer interclump gas could then possibly have 
appreciable emission compared to that of the clumps and reduce Xf appreciably. However, 
if there is true clump/interclump pressure equilibrium and if the density of the interclump 
gas is lower than that of the clump gas by an order of magnitude, then the temperature 
would be higher by an order of magnitude. This would bring down the optical depth of 
the ^^CO J = 1 — > line in the interclump medium by two orders of magnitude and could 
weaken its emission below what it would have been if Tj^ (intcrclump) = Tj^ (clump). If so, 
then the interclump gas would have an even weaker effect on the X-factor than has been 
estimated above. 

Of course, the above analysis assumes that the dominant CO-emitting gas is in the 
clumps, so that the presence of the interclump gas only increases Xf. However, if the 
interclump gas were to dominate the mass of the cloud and its ^^CO J = 1 ^ emission, 
then the appropriate average density to use in the expression for Xf would be considerably 
lower, thereby lowering Xf. 
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In summary, interclump gas would likely increase the X-factor by less than about 40%, 
so long as this gas did not dominate the mass of the cloud. If the interclump gas did 
dominate the mass, then the X-factor would be decreased, because the average density of 
the dominate CO-emitting gas would decrease. One other important effect of the interclump 
medium would be to increase the effective e. Estimates of e from inferred Xf values (and 
assumed or estimated clump densities and temperatures, see Section [3.2.1l) or from ^^CO J = 
1 0/^^CO J = 1 — s> line ratios (see Sections 13.31 and l4.4p could be higher than those 
expected from just clumps alone. 



3.3. ^^CO J = 1 ^ O/i^co J = 1 ^ 

One implication of the proposed formulation is that there is an additional layer of 
complication in interpreting line ratios. For example, if two spectral lines have t^j <^ 1, 
then the ratio of their intensities depends on the ratio their respective i7i/(T^) A(ro) values 
(where T^, is the excitation temperature of the transition), instead of depending on the ratio 
of their Jy{T^) [1 — exp(— tq)] values. If the two spectral lines have similar optical depths 
(i.e. similar A(ro) values) and Tg/ ^ 1, then thi s complication is rainimi zed. Given that the 



^^CO to ^^CO abundance ratio is about 60 (e.g.. iLanger &: Penziaslll990l ). the ratios of ^^CO 
lines to those of ^^CO still have this complication. 

To appreciate this, we start with a uniform slab of gas in which the J = 1 — lines of 
i^CO and ^^CO are in LTE: 

/(^=^C0) 1 - exp{-Xr.Ti2) 

IC'CO) ^ 1 - exp{-m) ' ^ ^ 

where Xr is the ^^CO/^^CO abundance ratio, X(^^CO)/X(^^CO). ri2 is the optical depth of 
12C0 J = 1 ^ and, in LTE, the optical depth of ^^CO J = 1 0, ri3, is related to m by 
Ti3 = Xr-Ti2. Expressiou fl57j) is an approximation because the integrated intensities are used 
here instead of the radiation temperatures at a particular velocity within the line profile. (It 
is also an approximation because the temperature corrected for the cosmic background and 
for failure of the R-J approximation, J'^{Tj^), is slightly different for the the two J = 1 ^ 
lines.) This line ratio has the following limiting cases: 

jf^UQQ-^ - ' < 1 and ri3 < 1, (58) 

~ 1 - eXp{-XrTi2) , for Ti2 > 1, (59) 
~ XrTi2 , for Ti2 ^ 1 ttud Ti^ <C 1, (60) 

~ 1 , for Ti2 ^ 1 and ri3 ^ 1. (61) 
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Equation fl58l) represents the uncommon case of optically thin ^^CO J = 1 ^ and has the 
unsurprising result that the line ratio is the abundance ratio when both lines are optically 
thin. Equation fl^U]) shows a linear relationship between the line ratio and the abundance 
ratio. When both lines are optically thick, as represented by equation flM]) . the line ratio is 
completely insensitive to the abundance ratio. 



If we now consider a clumpy medium and use the MSH84 formalism, then 

/(i^CO) 1 - ea;p(-are/,i2) 



/(12C0) 1 - exp(-re^,i2) 

;„ j-i -cx„„i-; j-;„„T J ii, 12nr 



(62) 



where Tg / 12 is the effective optical depth of the ^^CO J = 1 ^ line and where. 



a ^ ^ = . (63) 

Te/,12 A(ro,i2) 

Here it is implicitly assumed that A{tq) has the same functional form for both ^^CO J = 1 ^ 
and ^^CO J = 1 0. This is not nece ssarily the case gi ven that ^^CO can be selectively 



photodissociated more easily than CO (IWarin et al.lll996l ). Nevertheless, for simplicity, the 



same A{tq) for both of the J = 1 — > lines will be assumed here. Even for the same A(ro), 
the kA and e values can be different if ^^CO J = 1 — and ^^CO J = 1 ^ are in different 
optical depth regimes (i.e. tq^ 1 versus ro^3): 



I^A,12 ^0,12 



"'AA-.i '0,13 /^.N 

« = T. 7^ ■ (64) 



When ^^CO J = 1 — > is optically thin in the clump (i.e. ro,i3 <ti 1), /i;a,i3 = 1 and 613 = 1. 
If, at the same time, ^^CO J = 1 ^ is optically thick in the clump (i.e. to,i2^3), then it 
is easy to see that 

a = Xr Tq"^ , (65) 

where the subscript "12" was omitted from the k^ and e. In the other limit, when ^^CO J = 
1 — s> is optically thick in the clump, kA,iz = ^^,12 and eis = eu and 

a = xl , (66) 

where e is both €13 and ei2. With these expressions in mind, the following limiting cases 
result: 



/("CO) 
J(i2C0) 



Xr , for Tef^u < 1, Tef^i3 < 1, ro,i2^ 1, and ro,i3^1, (67) 

Xr k^^ tI)^^2^^ , for Tef^u < 1, Tef^i3 < 1, ro,i2^ 3, and ro,i3< 1, (68) 

K , for re/,12 < 1, re/,13 < 1, and ro,i3>3, (69) 
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^ XrTQ^i2 £ ^ , for re/,12 > 1, T^f,vi < 1, and ro,i3<l, (70) 

~ </CA^0,12 ^ ' ^e/.12 > 1' < 1> O'^C? ro,13^ 3, (71) 

~ 1 , /or re/,12 > 1 and re/,13 > 1. (72) 

The fc^ and e are those for ro^3. Only in the last case is the line ratio completely inde- 
pendent of the abundance ratio. When ro,i3^3, J(^'^CO)/I(^^CO) oc which is the same 
proportionality that the X-factor has in the case when the X(CO) is not spatially varying 
(see first paragraph of Section [3.2.3p . This represents a very important contrast with the case 
of uniform gas: even when both ^^CO J = 1 ^ and ^^CO J = 1 —>■ are optically thick 
on the central sightlines of the clumps, the observed ^^CO J = 1 ^ 0/^^CO J = 1 ^ 
line ratio does not reach an asymptotic value of unity, unless r^f^is ^ 1. And when both 
Te/,12 ^ 1 and re/,13 ^ 1 for ro,i3^3, the line ratio saturates near an asymptotic value that 
depends on the abundance ratio and the clump fluffiness. 

Expression (1^ and its limiting forms, (1^ to (1721) inclusive, can be applied t o the 



^^CO J = 1 and ^^CO J = 1 ^ data of the Orion A and B molecular clouds (Wall 



20061 ). The line ratio /(^'^CO)/I(^^CO) is plotted against r^f^u in the panels of Figure [51 The 
'Tef,i2 was found from solving equation (j4]) for Tg/. The J'uiTj^) was determined from the 
obseryed 1,^ (140 Atm)/I,^ (240 /xm) ratio and the one- and two-component models described in 



Walll (I2OO6I ). The uncertainties in r^f, a^r^f), were estimated from 



a(r.,) = ^ . (73) 

where 77// = T^^/ J'i,{Tj^). Any attempt to model the points in Figure [5] must consider 
the large scatter. The points with re/^0.3 and /(^^CO)/I(^^CO)^ 0.3, for instance, have 
large vertical scatter due to their large vertical error bars (which are really twice as large 
as shown, see figure caption). Even if the vertical scatter for these points is not entirely 
real, the high I C^CO) /1{^^C0) ratios still seem to be likely — i.e., I C^CO) /iC^CO) ~ 0.6. 
This combination of high /(^'^CO)/I(^^CO) and low Te/ is difficult to explain, and is best 
represented by expression fl69|) . For Xr = 1/60, a line ratio of /(^^CO)/I(^^CO) = 0.6 requires 
e = 0.12, roughly consistent with clumps that are hard spheres. Hence these points could 
be explained by gas with a low filling factor (and, accordingly, low t^j) and with clumps 
approximating hard spheres. 

On the other hand, the vertical scatter is probably real in the quasi-horizontal band of 
points that extends from t^j ~ to Xg/ > 1 for /(^^CO)/I(^^CO)< 0.3 to 0.5. Accordingly, 
no single theoretical curve can adequately model this band; models that have two curves 
each are plotted in the panels of Figure [51 Before the two-curve models are discussed, the 
one-curve models are described. 
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3.3.1. Crude Modeling of I{^^CO) /1{}'^C0) versus r, 



e/,12 



The data points in Figure [5] are those positions in the Orion A and B molecular clouds for 
which the peak T^^CO J = 1 ^ 0) > 3 a, and 1^(140 /xm), 1^(240 /im), and 1(^3^0) are all 
> 5 cr for a total of 372 points (see IWallll2006l ). One set of re/,12 values were derived from one- 



component models, hereafter T^f i2, and another from two-component models, hereafter T^f i2- 
These models were necessary for est imating the source function, J'u(Tj^), from the 140 /im and 



240 /im DIRBE data (see details in IWallll2006l ). While both panels of Figure [5] show roughly 



the same basic trend of I {^^CO) / 1 {^^CO) with re/,12, the models fit the I {^^CO) / 1 (^^CO) 
versus T^fi2 plot much better (i.e., upper panel) than the /(^^CO)//(^^CO) versus T^^^g (i-^-, 
lower panel). This is due to the unrealistically small errors in T^f 12- -^^^ these errors are only 
the formal errors. While the two-component models for the Orion clouds are better than the 
one-component models in many ways, they do not seem to adequately estimate the source 
function, j7^(Tj^). This source function depends on a model parameter, cq, that is not well 
constrained (unless combined with another parameter). Therefore, when systematic effects 
are considered, the effective errors in T^f i2 would be much larger. This must be borne in 
mind when considering the models described below. 

One of the simplest models to consider is that where each model curve is characterized 
only by a value of a, where a is defined in expression fl63l) . This kind of model represents 
a picture in which the clump properties are constant (i.e. fixed tq, Avc, Nc) and r^f^u and 
Tgj 13 increase only because the average number of clumps per sightline within each clump 
velocity width is increasing (i.e., increasing ^■^), thereby increasing the beam-averaged 
column density A^. In the limit where and aTef^u 0, (!62|) becomes 

^ a + 0.5 a (1 - a) re/,12 (74) 

If a is constant, then the line ratio would have a linear dependence on Te/,12 with slope 
0.5a(l — a) and intercept a. It turns out that ( |7il) is very good approximation even when 
the limit re/,12 — > is not valid. For example, (1741) yields a result that is within 2% of that 
of expression fl62|) even for re/,12 = 3, when a = 0.3. A sample of clumpy gas with identical 
clumps, but varying re/,12 would have a clear signature in a diagram like Figure [5j a straight 
line extending to the right from the intercept on the vertical axis with a positive slope, and 
there would be a well-defined relationship between the slope and the intercept. Figure [5] 
shows that, for the Orion clouds at least, such a simple signature is not obvious. 

Nevertheless, m odels with one curve were attempted using the orthogonal regression 



method described in IWalll (120061 ). Models with two curves were also tried. Neither the one- 



curve nor the two-curve models were successful. The two-curve models produced fits with 
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reduced chi-square values of xl = 13.9 and 24.4 for the T^f]i2 and the Tef,i2y respectively. 
These are factors of roughly 3 better than for the corresponding one-curve models. The 
recovered a values for the T^fi2 are 0.13 and 0.24 from the two-curve models. These are the 
same for the T^fi2- 

Models that are slightly more complicated are those that assume constant A(ro,i2)/Tcf,i2. 
The model curves here are characterized by two parameters: the A(ro,i2)/Tcf,i2 ratio, b, and 
the clump fiuffiness, e. This kind of model represents a picture in which the average number 
of clumps per sightline within each clump velocity width is constant (i.e., constant ^■^), 
but ro,i2 and tq i^ increase, thereby increasing Tgj 12 and Tgj 13. Combining equation (fTOl) with 
the constant b yields, 

T-0.12 = b^^ k~//' T% , (75) 
where kA and e are those for ro,i2 ^ 3. Substituting ( l65l) into ( 1751) gives 

a = b^ri'J^' , (76) 

where 

br = Xr k~^^'' b^ . (77) 

Substituting ( 1761) into ( l62l) results in 



/(i^CO) 1 - exp{-byj' 



. (7J 



/(12C0) 1 - exp(-re/,i2) 

The above expression results in growing very rapidly for small e. This growth cannot 

continue unchecked and must slow when a reaches its saturation value, given by expres- 
sion (!66|) . Consequently, in this limit, 

/(i^CO) 1 - ea;p(-<re;,i2) 



/(12C0) 1 - exp{-T, 



(79) 



e/,12 J 



Each model curve combines (1781) and ( 1791) : the minimum of the two /(^'^CO)//(^^CO) values 
for each re/,12 is used. As a result, the /(^^CO)//(^^CO) versus Te/,12 curves are from ex- 
pression (1781) for lower t^j^u and expression (I79|) for higher t^j^u- This produces an abrupt 
transition at some intermediate r^f^u value (see curves in Figure [5]). This abrupt transition 
is not real and only appears in these model curves because the smooth transition of A{tq) 
from optically thin to optically thick behavior (i.e., for l^ro^3) has not been considered 
here. Nevertheless, these intermediate tq values represent only a small portion of the curves. 
Accordingly, the abrupt transitions in the curves only introduce small inaccuracies into the 
final results. 
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Models with one curve and models with two curves were fitted, each fit optimizing bj- 
and e. The best fits of the one-curve models had reduced chi-square values that were factors 
of 2.1 to 2.5 worse than those of the two-curve models. The data points in Figure O are 
not entirely independent, so the effective number of degrees of freedom is about a factor 



of 9 lower than the number of points (see discussion of this in IWalll l2006l ). or about 40. 
Therefore, the F-test states that the two-curve models are better than the one-curve models 
at a 97.5% confidence level. The resultant parameter and reduced chi-square values are given 
in Figure O Even if the two-curve models are superior to the one-curve models, the fits of 
the former to the data are not impressive. Nonetheless, the fit to the T^fi2 is not too bad 
[xl — 4) considering the vertical spread in the points, even in the quasi-horizontal band. 

(2) 

The fit to the T^f 12 is considerably worse, but then this is due to the unrealistically small 
horizontal error bars. These two-curve models (where each curve uses two parameters) are 
better than the previous simpler two-curve models (where each curve uses one parameter) at 
the 95% confidence level or better, according to the F-test. For Te/,12 beyond some threshold, 
the constant a models (with one parameter per curve) superficially resemble the constant 
A(To,i2)/Tef,i2 modcls (with two parameters per curve), because the latter have constant a 
(i.e., a = x,^) in this r^f^u range. Consequently, the latter being significantly better than the 
former is because the latter models have a rapid non-linear rise at smaller Tg/ 12 where the 
majority of the points lie. 

The systematic uncertainties in the resultant parameter values were found by adopting 
systematic errors of 30% for /(^'^CO)//(^^CO), rjjf, and Xr- Specifically, 0.7 and 1.3 were 
each multiplied by one of these quantities and the models were refit, a total of 6 refits. The 
minimum and maximum br and e were chosen from these refits as measures of the uncertain- 
ties. The uncertainty in br is about a factor of 2 and for e about 30%. Expression (1771) and 
e ~ 0.3 imply that the systematic uncertainty in b is roughly a factor of \/2 or about 40%. 
Given that 1/6 is in reality this latter also has a systematic uncertainty of about 40%. 

With these systematic uncertainties in mind, the e values given in Figure [5] for the 
upper panel are consistent with those in the lower panel. The b^ values given in Figure [5] are 
significantly different between the two panels. Nonetheless, the changes in e from one panel 
to the other conspire with the br values to produce b values that are much closer together 
(in a logarithmic sense). Solving fl77|) for b, adopting kA = 1-5, and using the br values in 
Figure [5] results in 

1 X 10-2 and 6 x lO'^ , for r^^^, (80) 
Nc Av ^^'^2 ^ ^ 

= 2 X 10-2 and 8 x 10"^ , for rf^^^. (81) 



Considering that the effective spatial resolution of the observations was slightly less than 
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1° (or a linear resolution of 8pc for the adopted distance of 450 pc), these numbers are 
reasonable. The fluffiness values are more succintly summarized as roughly the following: 
e ~ 0.3 ±0.1. This range is consistent with that found in Section [3.2.11 

One concern raised by the modeling here is that many positions seem to have non- 
negli g ible o ptical depths in the ^^CO J = 1 — line. The modeling of physical conditions by 



Walll (120061 ) required ^^CO J = 1 — to be optically thin. Examination of the upper panel 
of Figure \5\ suggests that the majority of points are not towards optically thick sightlines 
in ^^CO J = 1 — 0. Most of the points are associated with the curve with the lower br 
value and, hence, lower A(ro,i2) values. The majority of these points are at lower Te/^u than 
occurs for the abrupt transition at re/,12 = 0.34. It is easy to compute that tq^is ^ 2 at 
this transition. Accordingly, the majority of points have tq i^ < 2. Given that tq i2 oc Teji2 
(cf. equation [75]), the ro,i3 is very strong function of Te/,12 and drops rapidly with re/,12 
decreasing below the transition point. Therefore, the majority of points do not represent 
sightlines optically thick in ^■^CO J = 1 — ^ 0, although some important minority of points 
may indeed represent such sightlines. 

In summary, the modeling of the J(^^CO)/I(^^CO) versus Tg/ 12 diagram did not produce 
exceptionally good fits, but does strongly favor one scenario over another. Specifically, it 
favors the case of growing clump optical depth as the column density increases from position 
to position, rather than the case of a growing number of clumps per sightline. The model 
fits suggest that the average number of clumps per sightline per clump velocity width is a 
few X 10~^ and the fluffiness is 0.3 ± 0.1, consistent with that expected for the observed Xf 
value (see Section [3.2.11) . 



4. Discussion 

The Orion clouds have been used as a test bed for the ideas developed here. More 
general implications are examined in this section. 



4.1. Comparisons with Previous Work 



Sakamoto ( 


1996 


) and 


Weiss et al. 


( 


2001) 



relationship between Xj and physical parameters — i.e., n(H2) and Tj^. The former used 
LVG models to explore the behaviour the X-factor under non-LTE conditions. The latter 
used observational data to estimate the dependence of Xf on density and temperature. 
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Sakamotd (119961 ) has the advantage that it deals with non-LTE conditions, but has the 
disadvantage that it does not consider radiative transfer through a clumpy cloud. Obviously , 
the advantage and disadvantage of the current paper complement those of ISakamotd (Il996l ). 
Accordingly, a treatment that combines non-LTE conditions with radiative transfer in a 
clumpy medium would be desirable, but is beyond the scope of the current work. 

Any attempt to at least roughly estimate th e properties of su ch a combined treatment 
must compensate for some minor errors found in ISakamotd (119961 ) . 

• The expression for the virial mass, equation (5) of that paper, is missing a factor of |. 
This can be verified for the uniform density case by comparing with the derivation given 
in Section [22] of the current paper. This can also be verifi ed for the p{r) oc case by 
comparing with the expressions given in iMaloney Jl990h. Both of the se comparisons 



confirm the absence of the factor ^ in equation (5) of ISakamotd (Il996l ). So equations 



2 

using the parameter must use ^ place of k^. Consequently, the numerical values 
of Xf computed from equation (7) of that paper must be corrected downward by a 
factor of -\/2- 

The claim that Xf is independent of Tj^ in the optically thin (i.e. r(^^CO J = 1 ^ 0) ^ 
1), thermalized limit is incorrect. This claim and its explanation are given in point ^^2 
(which contains expression 11) of ISakamotd (119961 ). Section [3.1.2! shows that such a 
claim is false. A very simple physical argument makes it obvious that Xf must depend 
on Tj^ in the optically thin, LTE limit: as Tj^ increases in the high-Tj^, LTE limit the 
CO molecules populate more and more of the upper-J levels, thereby decreasing the 
column density of CO in the J = 1 level. If the J = 1 — line is optically thin, 
then I{CO) decreases for fixed total N{H2), resulti ng in in c reasin g with increasing 
Tj^ (cf. Section [3.1.2p . The panels of Figure 3 of ISakamotd (1l996l ) do indeed show 
that Xf varies by less than factor of about 2 (from smallest to largest) over the range 
Tj^ = 10 t o 100 K for CO a bundances less than about 10~^. However, contrary to 
point #2 of ISakamotd (Il996l ). even n(H2) as high as lO^cm"^ does not thermalize the 
^^CO J = 1 ^ line when it is optically thin — LTE still does not apply. Given that 
the critical density of the ^^CO J = 1 ^ line is ~ 3 x 10^ H2cm~^, it is natural to 
assume that densities corresponding to the lower two panels of Figure 3 of that paper 
— i.e., 3 X 10^ H2 cm~^ and 1 x 10^ H2 cm~^ — are sufficient to bring the J = 1 — ^ 
line of CO close to LTE. The problem is that critical densities are normally computed 
as though a given atom or molecule is just a two-level system. For the CO molecule, 
there are many rotational levels and the radiative decay rates increase very rapidly — 
roughly like — as one moves up the rotational ladder. When the CO J = 1 ^ line 
is optically thin, radiative trapping effects are minimized and the CO molecules pile 
up mostly in the J = and J = 1 states, often resulting in superthermal J = 1/ J = 
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population ratios (i.e., the excitation temperature of the J = 1 ^ transition is 
greater than Tj^) or even population inversions, because of the higher levels with large 
radiative decay rates. Therefore, non-LTE effects play a role in the behavior of the 
J = 1 ^ line even at densities around 3 x 10^ H2 cm~^ (provided that the J = 1 
li ne is optica l ly thi n). For the specific cases of the densities considered in Figure 3 
of ISakamotd (119961 ) and for CO abundances ^ 10~^, the increasing Tj^ populates the 
J = 1 level from the J = level at an increasing rate, but also (depopulates the J = 1 
level to the higher- J levels at an increasing and roughly comparable rate. Accordingly, 
the column density of CO molecules in the J = 1 level stays roughly constant. Given 
that the CO J = 1 ^ line is optically thin for these abundances, the observed line 
brightness is directly proportional to the column density in the J = 1 (and not the 
J = 0) state, resulting in roughly constant l i ne bri ghtness and, ergo, roughly constant 
Xf. Therefore, expression (11) of ISakamotd (Il996l ). which states that Xf oc X{CO)~^ 
with no dependence on Tj^ nor n(H2), is approximately correct. And again, this is only 
in this particular non-LTE case. In the LTE limit when CO J = 1 — >^ is optically 
thin, the dependence of Xf on Tj^ is as given in Section 13.1.21 

The statement that the observed spectral profiles of the CO J = 1 ^ line place an 
upper limit on its opaci ty is erroii e ous, a nd is at odds with the clumpy cloud picture 
presented in Figure 1 of ISakamotd (119961 ). That paper states that the ^^CO J = 1 — 
line cannot be very optically thick or its profiles would be flat-topped. For a given 
volume of gas with more or less uniform density, that would be true; increasing the 
optical depth would saturate the line core before it would do so in the line wings. 
But this is not necessarily the case in a clumpy cloud. The optical depth of the 
^^CO J = 1 ^ line within each clump can increase without bound and, as long 
as the filling factor within each velocity interval within the line profile is sufficiently 
small, the observed line profile can still be Gaussian (e.g., see equation [7] of the current 
paper). This point is important for avoiding unrealistically low upper limits on the 



optical depth of the ^"'CO J 



line. 



With these caveats in mind, a relatively simple treatment of clumpy clouds that includes both 
radiative transfer and non-LTE effects can be developed. This will be particularly useful for 
developing insights into the emission of high-dipole molecules like CS, whose transitions are 
farther from LTE than those of CO. 

In the meantime, a rough estimate o f the b ehavior of Xf in the non-LTE limit is possible 
from combining the results oflSakamoto (119961) with the current work. For example, in the 
high-density, metal-rich case ISakamotol Jl996l ) finds that Xf oc n(H2)^/^ '^k^' ■'■^ current 
paper, this can correspond to any of the cases in Sections 13.1.31 13.1.4[ or I3.1.5[ All of these 
cases show a weaker dependence of Xf on n(Il2) and on Tj^. The side-on Gaussian cylinder. 
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for instance, requires roughly the following: Xf oc Cr T~°'^ n(H2)°'^. Therefore, as has been 
noted previously, the radiative transfer through a clumpy medium has the effect of softening 
the dependence of the X-factor on density and kinetic temperature. Two other cases are 
discussed in that paper. The case of intermediate density (i.e. ~ 3 x 10^ H2 cm~^) and low 
metallicity shows no dependence of Xf on n(H2) nor Tj^. Including radiative transfer through 
a clumpy medium is unlikely to change this. The case of low density (i.e. ^ 3 x 10^ H2 cm~^) 
and low metallicity has the following behavior: Xf oc n(H2)^^ '^k^^^ X(C0)~^. Combining 
this with radiative transfer in clump medium will yield Xf oc n(H2)^''^ 'j-<S2 X(CO)^^, where 
61 < 1 and 62 < |. Very crude guesses would be that 61 ~ 0.6 and 62 ~ 0.4. Something else 
worth considering is that in low-metallicity clouds, such as those in irregular galaxies, the 
CO abundance within the CO- emitting region within each cloud is close that in our Galaxy 
(see Ilsraellll997al . I2OOOL and references therein). Consequently, the behavior of the X-factor in 



such clouds may be closer to that described in Sections l3.1l and [3 . 2 . 2 1 than assuming optically 
thin, non-LTE ^^CO J = 1 — emission. 



Weiss et al.l (I2OOII ) observationally tested the dependence of the X-factor on density 



and temperature. They observed the central region of the galaxy M 82 in the lower- J lines 
of ^^CO, ^'^CO, and C^^O to constrain estimates of the molecular gas density and temper- 
ature. Using their estimates of N(H2), n(H2), and Tj^, they plotted N(H2)/I(C0) versus 
n(H2)^'^^ T~^ for the observed positions and found a roughly linear relationship (see their 
Figure 11). Determining the densities and kinetic ternperatu res from line ratios with suffi- 
cient accuracy is fraught with difficulties; IWeiss et al.l (120011 ) do not include the error bars 
in this plot and it is easy to show that the total horizontal uncertainty of each plotted point 
(i.e. from smallest to largest value for each point) is roughly 1/3 the horizontal extent of 
the entire figure. Accordingly, a large range of dependences on temperature and density are 
not ruled out. In general, such observational tests of the behavior of the X-factor are very 
uncertain. 



4.2. The Dependence of Xf on CO Abundance 



The molecular gas in low-metallicity irregular galaxies has very low CO abundance over 
large volumes and has ro ughly the Galactic abundance level in small internal regions within 



the molecular clouds (see Ilsraellll997al . l2000l . and references therein). Based on this picture. 



Section 13.2.21 discussed the dependence of the X-factor on the CO abunda nce and found 
that Xf oc X{CO)^^, where 9 can be anywhere from 0.5 to 1.0. Ilsraell (I2OOOI ) examined the 
dependence of Xf on the metallicity in a sample of Magellanic irregular galaxies and had the 
result that Xf oc X(0)~^'^. If the CO abundance is largely determined by the O abundance. 
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then Xf oc X(CO)~^'^. If, however, the CO abundance is largely determined by the C 
abundance and, given that X{C) oc X(O)^ '' (llsraell I2OOOI . and references therein), Xf oc 
X{CO)~^'^. Instead, the CO abundance might be strongly affected by the photodissociation 
rate and would be inversely pr oportional to that rate. Given that the photodissociation 
rate is proportional to X{0)~^ (llsraell Il997al . and references therein), this results in Xf oc 
X{CO)~'^'^. Hence, there are three possibilities: 



Xf oc X{CO) 



-2.5 



(82) 



or 

Xf oc X{CO)-^-'^ , (83) 

or 

Xf oc X(CO)-°-* . (84) 

Expressions (l82l) and (183|) represent dependences that are steeper than one would expect even 
in the optically thin case (i.e., Xf oc X(CO)~^ °). Having such a steep dependence would 
require that other physical parameters like density and kinetic temperature also depend on 
the CO abundance. When the volume-averaged X{CO) is small, the interstellar radiation 
field that reaches the CO-containing gas is strong, resulting in high gas temperatures. This 
would produce a smaller X-factor at small X{CO) (based on the expressions of Section ISTTl) 
and the dependence of Xf on X{CO) would be less steep. On the other hand, the strong 
interstellar radiation field would preferentially photodissociate the low-density molecular gas, 
leaving the high-density gas. The low X(CO) would be associated with higher densities and 
higher densities mean higher Xf (see Section IXTl) . Higher Xf at low X{CO) means a steeper 
dependence of Xf on X{CO). Consequently, the effects of the higher temperature would 
compete against the effects of higher density, making it uncertain if such steep relations are 
possible or even likely. The most likely dependence is that of expression fl84p . And this is 
consistent with the range found in Section 13.2.21 



4.3. The Relationship between Xf and ^^^CO J = 1 ^ O/^^CO J = 1 







That a relationship exists between the X-factor and the ^^CO J = 1 — > 0/^^CO J = 1 
line ratio is n ot new and has been u sed to determine whe t her th e X-factor varies spatially. 
For example, iRickard fc Blita (119851 ) and jPaglione et al.l (I2OOII ) observed the J = 1 — 
lines of ^^CO and ^^CO in external galaxies and found that the /(^^CO)/I(^^CO) ratio was 
significantly different in the nuclear regions from that in the disks. This is evidence that 
the X-factor in galactic nuclei is also different from that in the disks. Based on the current 
work, how do changes in J(^^CO)/I(^^CO) relate to changes in Xfl This can be assessed by 
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using an equation analogous to flMl) for ^^CO. Equation flM|) assumes that ^^CO J = 1 ^ 
is optically thin (i.e., to,i2^1). Since ^'^CO J = 1 ^ is likely to be optically thin on the 
scales of many parsecs, the appropriate expression is 

^f=CTK'^l~'^r'^-^^^ . (85) 

The Xf and K are those for ^^CO and not for ^^CO; accordingly the I{^^GO)/ I{^^GO) and 
are necessary to correct for those, respectively. Notice that when ^^CO J = 1 ^ 
is optically thin as well (to,i2^1), the line ratio is Xr (see equation WJ\} and equation ([H5D 
reduces to as required. When ro,i2 > 3, equation fl68l) for /(^^CO)//(^^CO) applies; 

when this is substituted into expression fl85|) . equation fl28|) for Xf is recovered, also as 
required. Therefore, the X- factor is proportional to the ^^CO J = 1 ^ 0/^^CO J = 1 ^ 
line ratio. 

However, the situation is not quite that simple. Given that the X-factor is also pro- 
portional to T^~^ x~^, determining spatial variations in Xf from one observed position to 
another can be obscured by unknown spatial variations in T^^ and x^- For example, Xf may 
be varying spatially while J(^^CO)//(^^CO) remains constant and vice versa. One way to 
estimate the dependence between Xf and /(^'^CO)//(^^CO) when there are also spatial vari- 
ations in Tj^ is to assume some kind of correlation between the Tj^ and n(H2). For example, 
Tj^ oc n(H2)"'. One of the limiting forms of equation fl62|) in Section [373] — i.e., equations fl67j) 
to flTTl) — is useful here: the most appropriate of these for the bulk of molecular clouds is 
equation fIBSl) . Using that expression along with that of (1221) and substituting into 
results in the following possible dependence: 



where 



/(12C0) 



(86) 



i(l-e) + wine -1.2,2) 



" = VTT ^ = 10 to 20 K. (87) 

(2 -7W^)(l-e) 

In the high-Tj^ limit, the "1.32" is replaced with "1.00". Of particular interest is pressure 
equilibrium, where w = —1. For e = 0.2 to 0.4 (see Section [3.3.ip . u = 0.8 to 0.7 when 
Tj^ = 10 to 20 K. In the high-Tj^ limit, the corresponding range in u is 0.5 to 0.3. Hence, 
if there are spatial variations of the molecular gas kinetic temperature from one observed 
position to another, and pressure equilibrium applies, the apparent correlation between the 
X-factor and the ^^CO J = 1 ^ O/^^CO J = 1 — * line ratio could be weak. 

There are, of course, additional complications. One is that equation fl85l) assumes LTE; 
non-LTE effects might be important in some sources and this would require an explicit density 
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dependence in the relationship between Xf and /(^^CO)//(^^CO). Another comphcation is 
that the ^^CO J = 1 ^ 0/^^CO J = 1 ^ ratio may be saturated (see equation [691) toward 
some positions in some sources. (See Section 1^741 for more details.) 

Nonetheless, if a large number of positions are sampled in a source where the kinetic 
temperature is roughly constant, such as in a galactic disk on large scales, then there could be 
an approximately linear correlation between the X-f actor and the ^^CO J = 1 ^ 0/^^CO J = 
1-^0 line ratio. 



4.4. Determining e from ^^^CO J = 1 ^ O/^^CO J = 1 ^ 

One way to constrain estimates of Xj is to model the data in the J(^^CO)//(^^CO) 
versus plot, as was done in Section 13.3.11 thereby estimating e and possibly other 

parameters. However, this method has th e considerable drawback that it requires some 



measure of the source function. IWalll (120061 ) did this for the Orion clouds by using the far-IR 
continuum data. The Orion clouds are 10° to 20° out of the Galactic plane, meaning that the 
origin of the continuum emission was unambiguous. Due to the lack of velocity information, 
continuum data is not very useful for isolating the emission of individual clouds along lines 
of sight through the Galactic plane. Therefore, determination of the source function of such 
clouds with continuum data would be very difficult, if possible at all. 

Another simpler approach is to use the histogram of /(^^CO)//(^^CO) values. One such 
histogram for the Orion clouds is shown in Figure [61 It shows a prominent peak between 
^^CO J = 1 ^ 0/^^CO J = 1 — s> line ratios of about 0.1 and 0.3 with a maximum at about 

0. 2. The left edge of this peak roughly corresponds to the self- shielding limit for molecules — 

1. e.. Ay ^ Imag. The ^'^CO J = 1 ^ 0/^^CO J = 1 ^ line ratio approximately correlates 
with the gas column density. For the Orion clouds, /(^'^CO)//(^^CO) = 0.1 corresponds to 
total hydrogen-nuclei column density, N(H I+2H2), of 0.7 to 1.5x 10^^ cm'^ or about = 0.5 
to 1.0 magnitudes. (Note that this range applies regardless of whether the column densities 
were derived from one- or t wo-compone nt models and regardless of imposed restrictions on 



the signal-to-noise ratio; cf. IWalll (120061 ).) The right edge of the peak can be interpreted by 
inspecting the panels of Figure 13 the value of J(^^CO)/J(^^CO) corresponding to the right 
edge of the histogram peak corresponds to where the model curves are linear — i.e., where 
the ^^CO J = 1 0/^^CO J = 1 ratio saturates. According to equation (1^ . this is x^. 
If Te is the ^^CO J = 1 ^ 0/^^CO J = 1 ^ ratio at the right edge of the peak, then e is 
given by 

_ ln{re) 
ln{xr) 
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If the ^^CO/^^CO abundance ratio, Xr-, is the same in the Orion clouds as in the inner 
Galactic disk, then the fluffiness, e, of the molecular cloud clumps are also the same. For 
Xr = 1/60, e = 0.29. This is cori s istent with model results of Section [3.3. 1[ However, the 



observations of iLanger fc PenziasI (Il990l ) show that Xr in the inner galaxy systematic varies 
with Galactocentric radius by a factor of about 2 from the center to the solar circle. Because 
of the logarithmic dependence on Xr the affect on e is relatively small: e = 0.29 to 0.35 from 
the solar circle to the inner Galaxy. And these values are still consistent with the model 
results of Section 13.3.11 even though those results only apply to the Orion clouds. Whether 
the X-factor is the same in the Orion clouds as in the inner Galaxy also depends on the 
typical densities and temperatures in the molecular cloud clumps in the two. If they are 
the same on average, then Figure H] provides a crude estimate of the X-factor: Xf ~ 1.2 to 
2.2X20 depending on the curve and e value (i.e. 0.29 or 0.35) used. 

Of course, not all positions have ^'^CO J = 1 — 0/^^CO J = 1 ^ ratios between 
0.1 and 0.3. A few are much larger, up to 0.95. Line ratios this high can only occur if 
either e is smaller than the value (s) that occurs for the positions in the peak, as discussed 
in Section 13.31 or if Te/^is is large. It is possible that there is a distribution of e values in 
clouds. This distribution would probably have to be sharply peaked to permit the peak in 
the histogram to have reasonably defined right edge (e.g., the peak from its maximum down 
to some poorly defined zero level on the right side has ^^CO/^^CO line ratios from about 0.2 
to about 0.45). To have a rough guess at what such a distribution might look like, it could 
be assumed that all the observed positions have saturated ^^CO J = 1 ^ 0/^^CO J = 1 ^ 
ratios. Then equation flSSl) with all the observed line ratios used in place of r^. will give a 
histogram that would contain such a distribution within its boundaries. Adopting symmetry 
about the maximum of the e distribution as a working assumption implies that there is an 
excess in the number of positions with high e values (e > 0.3), which corresponds to an excess 
in the number of positions with low ^^CO J = 1 — 0/^^CO J = 1 — ratios (< 0.3). This 
excess supports the interpretation that these low line ratios are unsaturated as expected, 
and do not represent part of this hypothetical distribution of e values. 

Comprehending exactly what limits the range of e values in a molecular cloud, cloud 
complex, or galaxy is key to understanding the X-factor. It is also key to comprehending 
the narrow range of ^^CO J = 1 — > 0/^^CO J = 1 — > ratios. Therefore, understanding the 
narrow range of ^^CO J = 1 ^ 0/^^CO J = 1 ratios is key to understanding the X-factor. 
These two are linked. Even if the formulation for the X-factor proposed in this paper were 
incorrect, even if it were not possible to clearly define a quantity like e in real molecular 
clouds, the connection between the restricted range of ^^CO J = 1 0/^^CO J = 1 
ratios and the X-factor would exist independently. 
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4.5. Clump/Interclump 



Real molecular clouds (like Orion and Rosette) have about | of their mass and of their 



volume in the form of clurnps, and the remain der in an interclump medium (see, e.g.,|Ostriker 



1999; Williams et al. 



the latter (as done in 



995 



Ba 



Ostrikei 



V et al.lll987D . Adopting 6h as the former quantity and fn as 



19991 ). then tm = 



H 



1 



I and To 



[(1 - fH)/fH] 



For Tm < 1, we need r„ ^ 



for the large limit to apply; since 2 ^ 0.7 is more or 



less valid, Xf would be increased by about (1 



'^mT'^^^~'^\ which is an increase of less 
than 15%. Therefore, if most molecular clouds are like the examples given here, then the 
interclump medium will have a small effe ct on the X-factor. This is especially true if this 
interclump gas is mostly atomic (e.g., see lWilliams et al.lll995l . and references therein). 



4.6. Radial Density Profiles 

The density dependence on radius in dark cl ouds or in the clurnps (or cores) of clouds 



is often p oc r where a is 2 or nearly so (e.g., \ 


iarvev et al. 


2001; 


Tachihara et al. 


2000; 


Lada et al. 


1999; 


Alves et al. 


1998; 


Henriksen et al 


|l997; 


Williams et al. 


1995). As discussed 



in Appendices O and iDl ( see Sections IC. 41 and [P. 3. 1|) and Section [3. 1.51 a p oc dependence 
requires specifying an outer radius, ri, to keep the total mass finite and an inner radius to the 
region, ro, to keep the density finite. Also, the ratio ti/tq determines the fluffiness, e, of 
the clump. And estimates of e can be constrained by observations Xf or of the ^^CO J = 1 ^ 
0/^^CO J = 1 — > intensity ratio (cf. Sections 13.31 and fS.S.ip . Accordingly, the observations 
can constrain the ti/tq ratio (assuming that the bulk of the ^^CO J = 1 ^ emission 
originates in structures with p oc r^^). The histogram of J(^^CO J = 1 ^ 0)/I(^^CO J = 
1^0) values suggests that e ^ 0.3, the same as that found in the modeling in Section [3.3.11 
for the Orion clouds. If the range of e values found from that modeling are relevant to much 
of the Galactic disk, then e = 0.2 to 0.4. If the clumps are spherical, then ti/tq would be 
in the range 2 to 9. If the clumps are cylindrical filaments viewed side-on, then ti/tq would 
be from 4 to 42. (End-on cylindrical filaments would have a much more restricted range of 
ri/ro, but these need not be considered here. It is much more likely to view high-aspect 
cyhnders close to side-on than close to end-on. Hence, the side-on cylinders would probably 
dominate the emission from most clouds.) 

Observations suggest ri/ro ratios that are consiste nt with the e values specified above. 



Specifically, ri/rn 


is found to be from about 4 to 15 (e.g.. 


Harvev et al. 


2001; 


Tachi 


lara et al. 


2000 




Lada et al. 


1999: 


Alves et al. 


1998: 


Henriksen et al. 


1997 




Williams et al. 


1995). It 



must be borne in mind that many of these observed Ti/tq values are merely lower limits due 
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to limitations in spatial resolution, which place upper limits on tq — the inner radius of the 
p oc region within each clump. 



4.7. Molecular Clouds in Some Specific Regions 



4^.7.1. Rosette and Orion 

The literature has papers with lists of clu mp temperatures, mas ses, and dimensions for 
certain molecular clouds. Two exa mples are IWillianis et al.l (119951 ) for the clumps in the 
Rosette Molecular Cloud (RMC) and lNagahama et al. I Jl998h for the clumps in the Orion A 
cloud. With these clump properties, rough estimates of Xf can be made. 

1995h 



For the RMC, Table 2 of Wilhams et al 



adjusting the adopted N {H2) / N {^^CO) ratio. IWiUiams et all (119951 ) adopt N(H',)/N(^^CO) = 
5 X 10 ^, which may be quite re asonable, especially fo r the pOph cloud (see iBertoldi fc McKee 
1992, and references therein) . Frerking et al. ( 19821 ) find an abundance of 8 x 10~^ for ^^CO 
and iLanger fc PenziasI Jl99oh find 1/60 for X(i3CO)/X(i2CO) in the Galactic disk approx- 
imately at the solar circle (and out to Orion). These imply N{H2)/N{^^C0) ~ 7.8 x 10^. 
Adopting this abundance ratio implies 55% higher clump masses (see M^te in their Table 2) 
and an observationally determined Xj that is also 55% higher or 1.7 X20, which is about the 
value expected for molecular clouds in the Galactic disk. At the very least, this higher 
N {H2) / N {}^C0) ratio is the same as that adopted for the Orion clouds. 



For the RMC clumps. IWilliams et al.l (119951 ) find 7^ = 5 (cf. their Figure 22), which may 



only be a lower limit due to resolution constraints. If taken at face value, this means that 
e = 0.30 and A;^ = 1.80. Section [Q4l of Appendix [Cl computes the specific expression for Xf 
for those parameter values — expression (1C43I) . The densities of the clumps were determined 
by assuming spherical clum ps with the masses (i.e ., M^te increased by 55%) and radii (i.e.. 



AR) as listed in Table 2 of IWilliams et al.l (119951 ). The densities and temperatures of the 
clumps either inferred from, or listed in. Table 2 of that paper will then give the Xf value 
corresponding to each clump — that is, the Xf value that would result in a hypothetical 
cloud filled with clumps identical to this one. The median of these Xf values is then taken 
as an estimate of the Xf of the cloud. (The mean is not used because it is influenced by 
extreme outliers. A better way to do this is to develop a formulation that includes clumps 
with a distribution of densities and temperatures, but that is beyond the scope of the current 
paper.) Here the Xf values use the density as averaged over the entire clump volume (i.e. 
equation IC46I1 and is more appropriate in this case (see IWilliams et al.lll994l . Il995l ). The 
Xf values found for the RMC range from about 2 to IOX20 with a median of 3.4X20- This 
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median is double that determined observationally (and adjusted to the new N{H2)/N{^^C0) 
ratio): 1.7X20- This suggests that the formulation is moderately successful, for the RMC at 
least, because it estimates the X-factor to within a factor of 2. 

It may be possible, however, to improve the agreement betwe en the model Xf and 
that observed. At least some of the temperatures listed in Table 2 of IWilliams et al.l (119951 ) 
are probably underestimated. These temperatures are excitation temperatures of the J = 
1 — s> transition of ^^CO that were computed from the peak radiation temperature, Tj^, of 
that line. In practical terms, this excitation temperature is about the sam e as the kinetic 
tempe rature, given that the ^^CO J = 1 — > line is easily thermalized (and IWilliams et al. 



( 119951 ) themselves implicitly assumed equivalence between the two temperatures when they 
determined column densities of all the molecular gas and not just the column densities in 
the J = 1 state of ^'^CO.) This method of determining excitation temperatures or, roughly 
equivalently in this case, kinetic temperatures implicitly assumes that the emitting gas fills 
the beam within each velocity interval of the line profile, especially at the line peak. Despite 
the scatter, the inferred clump temperatures have a definite positive correlation with the 
clump diameters; the Spearman rank-order correlation test places a 99.998% confidence level 
on the correlation (or the significance of the null proposition of no correlation is ~ 2 x 10~^). 
In contrast, in the Orion A cloud, where the kinetic (excitation) temperatures were als o 
estimated using the peak Tj^ of the ^^CO J = 1 ^ line (see iNagahama et al.l Il998l ). 



the Spearman test places only a 40% to 60% confidence on the correlation (depending on 

whether the clump size used was filament length or filament diameter). Accordingly, any 

correlation between clump size and inferred clump temperature in Orion A is either weak 

or non-existent. (Although it should be mentioned that these temperatures are not entirely 

reliable either, given that certain values are frequently repeated.) The Orion clumps have 

linear sizes (i.e., in parsecs) factors of 3 to 5 larger than those of the RMC and the Orion 

clouds are nearly a factor of 4 closer than the RMC. The interpretation of why the RMC 

clumps clearly show a correlation between size and inferred temperature while the Orion 

clumps do not is then fairly straightforward: the smaller clumps in the RMC are not being 

resolved in the velocity interval at the line peak, while the larger clumps are better resolved. 

Hence it is likely that at least some of the inferred clump temperatures for the RMC are too 

low. These inferred temperatures range from about 5 to 31 K with a median of about 9K. 

For Orion A, the range is 13 to 37 K w ith a median of about 18 K, w hich is the dominant dust 

temperature of the Orion clouds f see lWall et al.lll996l : IWallll2006l ) and even of the Galactic 

plane clouds (jSodroski et al.lll994l ). Both Orion A and the RMC are GMCs and should have 

clumps with roughly similar properties; the median temperature of the clumps of the RMC 

should not be half of that of the Orion A clumps. In short, there are three reasons why some 

of the clump temperatures in the RMC are probably underestimated: 

1. The temperature versus clump size correlation argues for a spatial resolution effect. 
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2. Both the RMC and Orion A are GMCs forming massive stars and their clumps should 
have similar temperatures. 

3. An additional reason is that the RMC is close to the Gala ctic plane (i.e. about 2°) and 



the dust temperatures in the plane is about 18 K (e.g., see lSodroski et al.lll994j ). If the 



KMC is like the Orion clouds in that its dust-gas temperature difference is small (jWall 



20061 ). then the gas kinetic temperatures must be roughly double the listed values for 
some of the clumps. 

Correcting all the RMC clump temperatures upwards by a factor of 2 would not be appropri- 
ate; the temperature-size correlation suggests that the correction factor should be largest for 
the "coldest," smallest clumps and progressively decrease towards unity for the "warmer," 
larger clumps. 

How do these corrections to the temperatures affect the inferred X-factor values? A 
simple way to address this is to raise all the listed temperatures by a factor of 1.5. This 
approach is less extreme than doubling all of the temperatures and is simpler than finding 
some prescription for applying different scale factors to the temperatures of different-sized 
clumps. This temperature correction will then affect the observed column densities and, in 
turn, the inferred averaged densities by a factor 1.5'^"^'^^ or 1.19. (Again, note that the 
"1.32" in the exponent becomes "1.00" for T,^ ^ 20 K.) The theoretical Xf will then change 
by a factor of 1.5"°-^^ x 1.19°-^^ or 0.88. The median Xf becomes 2.8 X2o- The observationally 
inferred Xf, however, increases by a factor 1.19 to about 2.0X20- Consequently, there is still 
a discrepancy between the Xf from the currently proposed formulation and that that would 
be inferred observationally; the theoretical Xf must be corrected downwards by 30%. 

There is a similar but somewhat smaller discrepancy for the Orion clouds. The Orion 



clumps are filaments (INagahama et al.l Il998l ) and are treated here as p oc r ^ cylinders. 



This type of cylinder was treated in Section ID. 3. II of Appendix [Dl but it was collapsing 
and magnetized. Here we assume the cylinders to be virialized and again adopt the aver- 
age density to be over the entire volume. This gives ky = 2.76 x 10~^^ in cgs units and 
kN = {''^i/ro)[2 — (ro/ri)][l + 2ln{ri/rQ)]~^. Here ti/tq = 14 is adopted because it has the 
fluffiness of that obtained from the modeling in Section 13 . 3 . 1 1 (also see Section numerical 
integration of ([2]) and comparison with (fTOj) results in e = 0.30 and kA = 1-76. Therefore, 

X;(X2o)=0.57CTT-°-^«nf^ (89) 



The Xf value corresponding to each clump listed in Table 2 of iNagahama et al.l (119981 ) was 
determined from the temperatures listed in, and the densities inferred from, that table. 
The densities were estimated from the listed filament masses and dimensions. For the Orion 
clouds, assuming cylindrical clumps and that the inferred densities are over the entire volume 
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will result in Xj values that range from 1.5 to 2. 7X20 with a median value of 1.9 X20, which 
is close to the observed value of ab o ut 2. 1 X^n ( lWallll2006l ). (Note that this Xf applies to 
the one-component models of [Wall (120061 ). The two-component models would imply that 



Xf ~ 3.5X20, ex cept that it is uncertain exactly how much mass is in the other colder 



component — see IWalll (j2006bl ).) One additional consideration is that at least some of the 



observed filaments are being viewed partly end-on. If we were to simplistically assume that 
there were only end-on filaments with the same dimensions, densities, and temperatures as 
those listed in Table 2 of lNagahama et al.l (Il998l ). then the median Xf would roughly double, 
even though the ti/tq = 14 cylinders considered here have e = 0.53 when viewed end-on. If 
we then consider that some of the filaments are partly end-on, then the best estimate of the 
theoretical Xf would be somewhere around 3X20- In fact, given that high-aspect filaments 
are more likely viewed side-on than end-on, this would probably be closer to about 2.5X20- 
Accordingly, the derived X-factor must be corrected downward by about 20-30% (assuming 
that ti/tq = 14 is appropriate for the Orion filaments). 

There are a few points to consider here. One is that average density defined as being over 
the whole cylindrical volume results in an X-factor that increases slowly with increasing e; 
this is quite different from the behavior seen with the previous definition of average density 
(see Section 13.2.11 and Figure H]). Another point is that defining the average density this 
way results in fewer of the Orion A clumps being virialized; the fraction of clumps with 
velocity widths within 40% of the virialized width drops from 74% to 51%. And the final 
point is that the derived X-factor values must be corrected downwards by about 20-30% to 
match those observed. This may only require a better definition for the average density. 
It is also possible that simply redefining the average density is not enough to account for 
the discrepancy between the theoretical and observed X-factor values. There must be some 
additional consideration to be included in the formulation. This is not surprising given 
that the current f ormulation simplistically neglects the effects of m agnetic fields and surface 
pressure (e.g., see lBertoldi fc McKeelll992l : iTilley fc Pudritzll2003l ). It is likely that improved 
understanding of molecular cloud physical conditions is also necessary for reducing this 
discrepancy. In any event, the derived X-factor is within a factor of 2 of that observed. 



4-7.2. The Galactic Center 



The molecular gas within the central few hundred parsecs of the Galaxy represents an 
environment that is distinct from that of the Galactic disk. Compared to the disk gas, 
the molecular gas in the Galactic cei iter is denser and hotter by an order of magnitude 



and very dynamically active (e.g., see iMartin et al.ll2004l ; iRodriguez- Fernandez et al.l 12001 
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1985 



Padione et all Il998l: iHiittemeister et all Il998l . Il993l : iBallv et all Il988l . Il987al : iHarris et al. 



Glisten et al. 19851 ). Consequently, it is no surprise that the X-factor for the Galactic 



center clouds differs greatly from that for the Galactic disk clouds. Specifically, it is between a 
factor of a few and more than an order of magiiitude smaller in the Galact ic center than in the 
disk (see, e.g.. lDahmen et al.lll998l : lOka et al.lll998l : ISodroski et al.lll995l ). S maller X-factors 



may be coramon in the central few hundred parsecs of spiral Ralaxie s fe.g.. iRickard fc Blitz 



1985 



Israellll988[ IWall et al.lEggsl : iReganlboOol : IPaglione et allboOlh (see also 



Dahmen et al 



19981 . and references therein). 



How do we account for these low X-factors? Given that tidal forces may be appreciable 
in the Galactic center, the outer layers of clump s or clouds may b e sheared off, producing 
a substantial interclump or intercloud medium (jStark et al.lll989l ). As mentioned in Sec- 
tion [2231 if interclump gas dominates the ^^CO J = 1 ^ emission and the gas mass, then 
Xf decreases because of the low average density in this gas. As stated in the previous para- 
graph, there is evidence for dense, hot molecular gas in the Galactic center region, but does 
this apply to the bulk of the CO-emitting gas in the Galactic center? If, for example, this CO- 
emitting gas is an order of magnitude hotter and an order of magnitude less dense than that 
in the Galactic disk, then the low Xj migh t be explained. However, the evidence for this is 
less than compelling. IPahmen et aP Jl998[ ) find high I{^^CO J = 1 ^ 0)/I(C^^O J = 1 ^ 0) 
ratios that are sometimes nearly as high as the ^^CO/C^^O abundance ratio. They argue 
that the higher the J(^^CO J = 1 -> 0)/I(C^^O J = 1 ^ 0) ratio, the lower the gas density, 
claiming densities as low as n(H2) ~ 10^ cm~^, even though a simpler explanation is that the 
higher ratio only really implies lower optical depths in the J = 1 ^ lines of ^^CO and ^^CO. 
These authors used observed molecular gas velocity gradients (or velocity widths per unit 
size) to impose constraints on the model velocity gradient. These observed gradients only 
really apply to the scale of entire clouds, and applying them to the model gradients is no 
more vahd than applying such constraints to gas volume density. Such cloud-scale quantities 
are only rough lower limits to the relevant quantities to be used in non-LTE radiative codes, 
such as the LVG code. Assuming that the lowest and highest observed cloud-scale velocity 
gradients represent both lower and upper limits on the velocity gradients to be used in the 
LVG code is not correct and imposes an artificial connection between optical depth and den- 
sity — both rise and fall together. Thus when the I{^^CO J = 1 ^ 0)/I(C^^O J = 1 ^ 0) 
ratio approaches the abundance ratio, the optical depths of both ^^CO J = 1 — » and 
C^^O J = 1 — s> drop and, because of the imposed constraint on the velocity gradient, the 
density drops too. Thus imposing constraints on the velocity gradients in this way leads to 
the assertion that I{^^CO J = 1 ^ 0)/1{C^^O J = 1 ^ 0) being close to X(i2C0)/X(CiS0) 
cannot occur in LTE — which is incorrect. Densities of n(H2) ~ lO^cm"^ may indeed be 
present in the CO-emitting gas in the Galactic centre region, but it cannot be verified with 
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this line ratio alone; additional information is needed (and something more than large-scale 
velocity gradients). 

Nonetheless, partly explaining the low X-factors in the Galactic center region is still 
possible. The large observed values of the I{^^CO J = 1 -> 0)/I(C^^O J = 1 0) ratio 
( jPahmen et al.lll998l ). for example, suggest that ^^CO J = 1 — is optically thin or nearly so. 
Prom Table [U we see that the optically thin case has Xf an order of magnitude or more lower 
than the Hard Sphere case, this latter having roughly the Xf for Galactic disk clouds when 
Tj^ = 20K and n(H2) = 2000 cm~^. One problem with this is that there is evidence that 
the molecular gas could have high kine t ic temperatures, e.g. T^^IOOK 



Harris et al 



2004J: iRodriguez-Fernandez et al.l l200ll : iHiittemeister et al.l Il998l . Il993l : iGiisten et al 



(e.g.. 


Martin et al. 


Glisten et al. 


1985; 



19851 



1^ If Tj^ is 100 K, then Xf rises to about 0.5 X20 (see equation 



m 



Section r3.1.2l) . or a factor of ~ 4 lower than that for the Galactic disk clouds. To have more 
than an order of magnitude lower Xf, an additional effect is necessary. The Galactic center's 
clouds are very dynamically a ctive and have line widths factors of 5 to 10 larger than found in 
Galactic disk clouds (see, e.g.. lBally et al.lll987al . Il988l ). If this also means that the velocities 
of the Galactic center gas are factors of 5 to 10 beyond virialization, then the effective of 
this gas is larger than the virialized value by the same factors. Given that Xf oc kl~^, Xf 
would be reduced by factors between 1 and 5 to 10, depending on e. If the ^^CO J = 1 — > 
line truly is optically thin, then e = 1 and the velocity structure is irrelevant (i.e., only 
iridirectly relevant because the velocity structure partly determines the optical depth) for 
determining Xf. Not all the observed /(^^CO J = 1 ^ 0)/I(C^®O J = 1 ^ 0) ratios are high 
enough to indicate optically thin ^^CO J = 1 ^ 0, so an intermediate e of 0.5 might be more 
appropriate. But this would only result in an X-factor that is only about a factor of 6 lower 
tha n the disk value. Fin ding X-factors that are more than an order of magnitude smaller 



see 



Dahmen et al.l 119981 ) in molecular gas that is warm and optically thin (or nearly so) is 
difficult to explain. Heightening the abundance of ^^CO by a factor of a few might do this, 
but it would lead to other difficulties, such as explaining how ^^CO J = 1 — could be so 
close to being optically thin (as suggested by the high /(^^CO J = 1 ^ 0)/I(C^^O J = 1 ^ 0) 
ratios). The best way to reduce Xf seems to be assuming that the gas behaves like fluffy 
clumps, something like the Squared Lorentzian Sphere. If the gas has a temperature of 
about 100 K and is still marginally optically thick, then it is easy to obtain Xf ~ 0.4 X20 — 
a factor of about 5 smaller than in the Galactic disk. If we also consider that the gas velocity 
widths are an order of magnitude larger than required by virialization, then this reduces the 
X-factor another factor of nearly 3 and it can be an order of magnitude smaller than that in 
the disk. 



In short, it is possible to explain X-factor values that are an order of magnitude smaller 
than is found in Galactic disk molecular clouds, but the really low values — nearly two 
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orders of magnitude smaller (see lDahmen et al.lll998l ) — are not so easily accounted for 



J^.1.3. Irregular Galaxies 



As discussed in Sections 13.2.21 and 14.21 the unusually high X-factors found in molecular 
clouds in irregular galaxies are l argely due to the low abundances of ^^CO (or low X(^^CO)) 
in these galaxies. As listed in [Israeli (119971 ). the Xf values found in molecular clouds of 
irregular galaxies can range from a factor of a few to about 100 higher than that found in 
the Galactic disk molecular clouds. As described in Section [3.2.21 the abundance X(^^CO) 
is an average over the volume of the molecular gas; there are CO-emitting regions with 
a Galactic abundance of CO surrounded by envelopes of CO-deficient molecular gas (see 
Israeli 1 199 71 . l2000l . and references therein). Consequently, a map of some molecular gas tracer 
that can show the molecular gas not traced by CO can be compared with CO maps to test 
the correlation betw een the observed Xf and the inferred volume ratio of molecular gas 
to CO-emitting gas. iMadden et al.l (Il997l ). for example, map the low-metallicity irregular 
galaxy IC 10 in the [C II] 158 /xm line and find that this emission is largely associated with 
molecular gas in photodissociation regions (PDRs). Their Figure 5 shows the contours of 
the [C II] 158 /xm-line emission with superposed contours of CO emission. Depending on 
exactly how the relative areas are estimated, the ratio of projected areas of molecular gas 
to CO-emitting gas is about 10 to 20. Assuming that the volume goes like (area)^'^, then 
these correspond to volume ratios of about 30 to 90. Accordingly, the volume- averaged 
X(^^CO) would be roughly the Galactic disk value divided by these volume ratios. Given 
that the X-factor goes like Xf oc X(^^CO)^°-^ (see S ection I4l2l). Xf would be factors of 
about 15 to 40 larger than found in the Galactic disk; iMadden et al.l (119971 ) find that this 
factor to be about 100. Therefore, using only very crude estimates of the volume ratio of 
molecular gas to CO-emitting gas, the X-factor is estimated to within a factor of a few. If 
the molecular gas in IC 10 is very roughly representative of that in other irregular galaxies, 
then it is easy to understand the high X-factors in these galaxies. If Xf/Xfc is the ratio 
of the observed X-factor to the "standard" value found in Galactic disk clouds, and if Ar 
is the ratio of projected areas of molecular gas to CO-emitting gas, then Xf/Xfc — Al^'^ is 
expected. Deviations from this dependence would be extremely interesting, providing clues 
to differences in physical conditions (e.g., e, n(H2), Tj^) between those in the molecular clouds 
of irregular galaxies and those in Galactic disk clouds. 
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5. Summciry and Conclusions 

An improved formulation for the X-factor is proposed that combines viriahzation of the 
gas with radiative transfer in a clumpy medium. The statement that the velocity-integrated 
radiation temperature of the ^^CO J = 1 — > line, /(^^CO), "counts" optically thick clumps 
is quantified using the formalism of MSH84 for line emission in a clumpy cloud. Adopting the 
simplifying assumptions of thermalized ^^CO J = 1 — > line emission and isothermal gas, an 
effective optical depth, Te/, is defined as the product of the clump filling factor within each 
velocity interval, and the clump effective optical depth (or effective optically thick 

area), A(to), as a function of the optical depth, tq, on the clump's central sightline. The 
A(ro) is well approximated (to within about 13-26%) as a power law in tq with power-law 
index, e, called the clump "fluffiness," and has values between zero and unity. While the 
^^CO J = 1 — s> line is optically thick within each clump (i.e., high tq), it is optically thin 
"to the clumps" (i.e., low Tg/). Thus the dependence of I{CO) on Tg/ is linear, resulting in 
an X-factor that depends only on the properties of the clumps rather than having a direct 
dependence on the entire cloud. Assuming viriahzation of the clumps yields an expression 
for the X-factor whose dependence on physical parameters like density and temperature are 
"softened" by power-law indices of less than unity that depend on the fluffiness parameter, e. 
The X-factor provides estimates of gas column density because each sightline within the beam 
has optically thin gas within certain narrow velocity ranges. Determining column density 
from the optically thin gas is straightforward and parameters like e then allow extrapolation 
of the column density of the optically thin gas to that of all the gas. Implicit in this 
formulation is the assumption that fluffiness is, on average, constant from one beam to the 
next. This is less true for density and temperature for which the X-factor, Xf, may have a 
weaker dependence. 

The proposed formulation addresses the problems of the explanation proposed by DSS86: 

1. Treatment of radiative transfer. The dependence on fluffiness, e, represents a radiative 
transfer parameter of the clumps. The optically thin case results by simply setting 
e — 1: dependence on viriahzation disappears and the expression for column density 
in terms of I {CO) in the optically thin case remains. In a hypothetical completely 
optically thick case with flat-topped line proflles, e — > and the radiative transfer 
vanishes, leaving only viriahzation (as in DSS86). 

2. Reduced sensitivity to and n(H2). The effect of the fluffiness reduces the dependence 
of Xf on density and temperature from Xf oc — to roughly Xf oc ^ "(j^^) j in 
the high-Tj^ limit. Thus variations of an order of magnitude in either n or Tj^ would 
allow X to vary by less than a factor of 2. 
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3. Virialization of entire clouds is unnecessary. The densities required to give reasonable 
values of Xf are consistent with those found in cloud clumps (i.e. ~ 10^ cm~^). 
Thus virialization of clumps, rather than of entire clouds, is consistent with the ob- 
served values of Xf. And even virialization of clumps is not strictly required; only a 
relationship between clump velocity width and column density similar to that of virial- 
ization can still yield reasonable values of the X-factor. The underlying physics is now 
at the level of cloud clumps, implying that the X-factor can probe sub-cloud structure. 

4. Stronger dependence of peakT^ on N(H2) than of Av on N(H2) is now explained. The 
peak Tj^ depends linearly on the filling factor within each velocity interval, 
thereby accounting for its dependence on N{H2). The Av, on the other hand, is the 
observed line width and is not necessarily directly related to the beam-averaged column 
density, N{H2). If virialization is important at the level of the clumps, then the velocity 
width-column density relation is between that of the clump velocity width, Avc, and 
clump column density, Nc. 

X-factor values were computed for both spherical clumps and cylindrical clumps (i.e. 
filaments) of densities 2 x 10^, 2 x 10'^, and 2 x 10^ H2 cmT^ and kinetic temperatures 10 and 
20 K and different internal density variations. The clumps of average density 2 x 10^ H2 cm~^ 
that reproduce the standard observed Xf value of about 2X20 for the Galactic disk clouds 
within a factor of 2 are the Hard Sphere (or uniform- density sphere), the Gaussian Sphere, 
and the Gaussian Filament (see Tabled]). Aside from the clumps listed in Tabled], spherical 
clumps and filaments with average densities of ~ 10^ H2 cm~^ and an density variation 
can also produce X-factors within a factor of 2 of the standard value. Testing these clump 
types reveals a potentially strong inverse dependence of Xf on e, as shown in Figure HJ 
Increasing the fluffiness has the advantage of weakening the dependence of the X-factor 
on density and temperature, but the disadvantage of decreasing Xf to values appreciably 
below the standard value. This strong dependence of Xf on e is related to how the average 
density is defined. For clump types with no clearly defined outer radius (e.g., Gaussian or 
Squared Lorentzian), the average density was defined as that within a spherical or cylindrical 
(depending on clump geometry) diameter equal to the FWHM of the projected surface 
density distribution. For clumps with an density variation, an outer radius must be 
defined to keep the clump mass finite. The average density can be defined as done previously, 
but also can be defined over the whole clump volume. The latter definition can result in a 
weakly rising Xj as a function of e. 

The proposed formulation also suggests a specific dependence of the X-factor, Xj, on the 
^^CO abundance, X(^^CO). In the molecular clouds in irregular galaxies, the CO abundance 
within the CO-emitting regions is about the same as that in Galact ic disk clouds, bu t the 



average over the entire molecular cloud volume is much lower (see llsraell Il997l . I2OOOI . and 
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references therein). If X(^^CO) is this volume-averaged abundance, then the formulation 
predicts Xf oc X(^^CO)'''^ *° '^ '^ (for e ^ 0.3). This is consistent with observations (see 
Section SSI). 

This formulation has implications for the interpretation of spectral line ratios, especially 
for the /(^^CO J = 1 ^ 0)/1{^'^CO J = 1 ^ 0) ratio. Modeling the plot of this ratio against 
the effective optical depth, r^f, — determined from the peak radiation temperature of the 
^^CO J = 1 — > line normalized to the source function (in temperature units) — for 
the Orion clouds provides crude estimates of the fluffiness: e ~ 0.3 ± 0.1. Histograms of 
^^CO J = 1 — * 0/^^CO J = 1 — *^ line strength ratios can also provide estimates of e, and 
have the very important advantage that no estimates of the source function are necessary. 
Such a histogram from observations of the Orion clouds shows a peak that extends from 
/(13C0 J = 1 ^ 0)/I(i2co J = 1 ^ 0) ~ 0.1 to 0.3. These limits for the Orion clouds 
(and possibly for other Galactic disk clouds) correspond to the minimum column density for 
self-shielding against the interstellar radiation field (i.e. ~ 1 mag) at the low end and to 
the saturation of the ^^CO J = 1 — 0/^^CO J = 1 ^ line ratio at the high end. The value 
of this ratio at saturation is determined from the dominant value of e (i.e. ~ 0.3) within the 
clouds' substructures. Consequently, a narrow range in e can simultaneously account for the 
limited range of ^^CO J = 1 — *• 0/^^CO J = 1 ^ line ratios and for a relatively constant 
X-factor. In any event, it is no surprise that ^'^CO J = 1 — »• 0/^^CO J = 1 — >^ line ratio is 
related to the X-factor. 

Observations of the ^'^CO J = 1 ^ 0/^^CO J = 1 ^ line ratio have been used, for 
example, to infer differences in the X-factor be tween the molecular clouds of the nucleus o f 



a spiral galaxy and that of its disk clouds (e.g., iRickard &: Blitzlll985l : iPaglione et al.ll200ll ). 
The formulation finds a linear relationship between Xf and /(^^CO J = 1 — > 0)/I(^^CO J = 
1 0), provided that Tj^ and X(^^CO)/X(^^CO) (or Xr) are constant. If these quantities 
vary spatially, then Xf can have a dependence on I(^^CO J = 1 — 0)/I(^^CO J = 1 — 0) 
that is weaker than linear. 

The proposed formulation and observed J(^^CO J = 1 ^ 0)/I(^^CO J = 1 ^ 0) ratios, 
or observed X-factor values, can constrain estimates of properties of substructures within 
molecular clouds, often by imposing limits on estimates of the fiuffiness, e. Many clumps or 
small clouds have an density (or p) dependence. If the clump outer radius is, ri, and the 
inner radius of the p oc region within the clump is tq, then the ratio ri/ro is constrained 
by the limits on e. Given that e ~ 0.2 to 0.4, and assuming that p oc structures are the 
dominant source of the ^^CO J = 1 — emission, spherical clumps would have ri/ro = 2 to 
9 and cylindrical clumps would have ri/ro = 4 to 4 2. Observations are apparently consistent 



with these limits, where ri/ro is about 4 to 15 (e.g.. lHarvey et al.ll2001uTachihara et al.ll2000 
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Lada et al.l Il999l : lAlves et all Il998l : iHenriksen et all 119971 : IWilliams et al.l Il995h . although 
spatial resolution limitations often mean that the observed numbers are merely lower limits. 

The properties of real clumps in real molecular clouds can be used to estimate the 
X-factor within these clouds and then be compared with the obs ervationally determined X- 
factor. Applyin g this to the Orion A c loud (INagahama et al.lll998l ) and the Rosette Molecular 
Cloud or RMC (IWilliams et al.lll995l ) yields estimates of the X-factor that are within a factor 
of 2 of the observed values. While this is acceptable as a start, reducing this discrepancy will 
require improving the formulation. Simply changing the definition of the average density can 
reduce the discrepancy, but it seems likely that something else is missing from the current 
formulation. 

A future, improved formulation for the X-factor must address the following shortcomings 
of the current formulation: 

A) It is not entirely clear why the fiuffiness seems to be constant, or at least sharply 
peaked, at one value. This value seems to be about e ~ 0.3. 

B) The resultant X-factor values are too dependent on the precise definition of the average 
density. A closer to optimal way of defining such a density must be found. 

C) The desired insensitivity of the X-factor to density and temperature often comes at the 
price of too low an X-factor value. Again, this is affected by the definition of average 
density. 

D) There is a potentially strong dependence on the fiuffiness, e. This dependence can be 
weakened by the appropriate choice for the definition of average clump density. 

E) The current formulation does not consider clumps with a spectrum of properties, such 
as distributions of densities and temperatures. This will introduce other parameters 
in addition to e, Tj^, and n(H2) for determining the X-factor. 

F) Additional physical effects must be considered. These would include non-LTE effects 
and the effects of magnetic fields, surface pressure, and turbulence. 

G) Observational determinations of parameters such as e can be complicated by clumps 
having a spectrum of properties, such as ranges in densities, temperatures, optical 
depths, etc. The presence of interclump gas might complicate this as well. 



Despite these shortcomings, the currently propose d formulation represen ts the first ma- 
jor improvement in understanding the X-factor since iDickman et al.l (119861 ) (i.e., DSS86), 
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because it includes radiative transfer. Previous explanations of the X-factor involved count- 
ing optically thick clumps (or entire clouds) and a relationship between the gas column 



density and the velocity width of the CO J = 1—^0 spectal line (e.g., DSS86 Ilsraellll988 



EvansI Il999l ). But applying the DSS86 approach directly to cloud clumps often overesti- 
mates the X-factor. At first glance, applying radiative transfer to an optically thick line is 
apparently pointless. However, portions of the gas are not optically thick, permitting esti- 
mates of the mass of the gas when radiative transfer is considered. And including radiative 
transfer does indeed result in reasonable estimates of the X-factor even when applying the 
formulation to the level of individual clumps. 

Determination of the X-factor on scales of many parsecs can constrain the average 
properties of the molecular gas at scales of just a few parsecs. Future formulations may 
refine the X-factor into a potent probe of molecular cloud structure. 
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A. The Effective Optical Depth of a Single Clump 

MSH84 derive the effective optical depth of a clump, A(ro), which they call the effective 
optically thick area, for the case of a clump with a Gaussian spatial variation of its optical 
depth. This appendix gives a quick derivation of A(ro) for a clump with a more general 
optical depth variation. MSH84 derive the effective optical depth on a line of sight through 
a clumpy cloud: 



'Tef{vz) = /^'^ ^ — / dv exp ( — ^ \ ^ ) I dx I dy <{! — exp 



2nAvJ " V 2 A 



-r(x, y) exp 



2Av2 

(Al) 

Notice that different notation from that of MSH84 is used here. The A/'c is the number 
of clumps per unit projected area of the cloud and corresponds to the N of MSH84. The 
velocity widths here are rms velocity widths, whereas those of MSH84 are ratios of the 
Gaussian line profile areas to their amplitudes. This results in extra factors of ^/2tt. The 
spectral line velocity width was represented by a in MSH84 and is given by A/27rAv here. 
The velocity width of a single clump was Vo in MSH84 and is -\/27rAvc here. The velocity 
Vz is the bulk velocity of a clump along the sightline and the velocity, v, is the velocity of a 
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given element of gas along the sightline. Following MSH84, the assumption of the line width 
being much greater than the velocity width of a single clump, i.e. Av ^ Avc, is adopted. 
This means that v cannot deviate too far from Vz without making the integrand very small. 
At the same time, the first exponential factor of the integrand does not change much from 
ea;p(— w^/(2 Av^)), because Av is so much larger than Avc. Consequently, expression (lAll) 
simplifies: 

T^,{vz) = -^J^ exp (-^) Jdv' Jdxjdyl^l- exp [-r(x, y) exp (-^) ] } 

(A2) 

where v' replaced v — Vz- Now the integrals are only over the clump velocity width and the 
clump's projected surface area. Therefore, the integrals, when properly normalized, must 
give the effective optical depth of a single clump, A(ro). Normalizing for the velocity width 
requires dividing by v^27rAvc. Normalizing for the effective projected area of the clump 
requires dividing by a^^^ (see equation [3]). This results in equation ([2]), as desired. 

Equation (IA2p then becomes 



a „A/'cAvc . _ / V 



2 



'./V"W = -^^AMexp^-^j . (A3) 

Given that a^^.^.Mc = NfNc, equation (1A3I) becomes equation ([T]), also as desired. 

Equation ([1]) can be used as the starting point instead. It represents a logical relation- 
ship between the effective optical depth through the cloud and that for a single clump; its 
derivation is trivial. Comparing equation ([T]) with equation (1A2I) then gives equation ([2]) for 
A(ro). 

The reader may have noticed that the definition of A(ro) is lacking a certain geometric 
correction. For example, in the tq ^ 1 limit, A{tq) ~ tq. And yet A(ro), which is an 
appropriately averaged r over the projected surface of the clump, cannot simply be the r 
through the center of the clump, i.e. tq. It must be less than this value, at least for those 
clumps where tq is highest on the central sightline. Consider, for example, the hard sphere, 
which means a sphere with clearly defined edges and uniform internal density. (See MSH84 
for their treatment of hard spheres.) Since we are only considering LTE and uniform kinetic 
temperatures, and that these clumps have uniform internal density, the optical depth on 
any given sightline through the clump is proportional to the path length through the clump 
on that sightline. The average path length through a hard sphere is its volume divided by 
its projected area or where d is the sphere's diameter. The average optical depth of 
such a sphere in the optically thin limit is then |ro rather than simply tq. However, the 
clump optical depth, Nc, as defined in the text is also used as the column density on the 
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central sightline through the clump. This means that N/Nc is not exactly the number of 
clumps in the beam as stated in the text: a slight correction factor is needed. Nevertheless, 
this correction factor cancels that needed for A(ro), thereby yielding the correct result for 
Tef. Another way to argue this is to consider the factors A(ro)/Nc and A^Avc/Av. The 
former factor is the average optical depth per unit column density for a clump, since both 
A(ro) and Nc are referred to the central sightline. The latter factor is the beam-averaged 
column density within velocity interval Av,, at the line central velocity, = 0. Therefore, 
the product of the the two factors must give the optical depth at the line central velocity, 
and this is equation ([1]) evaluated at = 0, as desired. 



B. The X-Factor for a Small Beam on a Uniform-Density Cloud 

Here the N{H2)/I{C0) conversion factor is estimated for the case of a beam much 
smaller than the source, where that source is a uniform- density cloud. The velocity-integrated 
radiation temperature of the ^^CO J = 1 — >^ line, I (CO), in LTE in a cloud of uniform 
kinetic temperature is 

I{CO) = MT^) [ dv [1 - exp(-r(v))] , (Bl) 

J line 

where the integral is over the spectral line profile. The optical depth velocity profile is 

r(^^) = roexp(-^) . (B2) 
Numerically integrating ( IBll) shows that the integral can be approximated by 

/ dv [1 - exp{-T{v))] ^ AvkAT^ , (B3) 

J line 

where = 1.3 and e = 0.14. If we now assume a uniform density sphere and virialization, 
then the derivation of Av of Section [2^2] is relevant. It is also assumed that the velocity width 
on whatever sightline through the cloud reflects the virialization of the entire cloud and does 
not have appreciable contributions from other types of motion. Repeating the derivation 
of Section [2.41 again yields equation fl28p . This means that this equation and its variants 
(see Appendices O and |D]) are more general than only treating radiative transfer through a 
clumpy medium. All that is required is that I{CO) oc Tq and Av oc (NL)°'^, where and L 
are the column density and pathlength, respectively, on the sightline through the cloud. In 
fact, the density and kinetic temperature that appear in expression fl28l) must be replaced, in 
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the current treatment only, by the corresponding quantities that are appropriately averaged 
over the entire cloud, rather than for just an individual clump. 



C. Treatment of Spherically Symmetric Clumps 

Here we assume that the mass density p is simply a function of the radius, r, within each 
spherical clump: p = p{r). This is similarly true for the number density, n(r). If the x-y 
plane is perpendicular to the observer's sightline and z is measured along the sightline, then 
n(r) is n(^/^~+~z^), where p = ^/x^ + y'^ is the projected radius as seen by the observer. 
The clump column density as a function of projected radius, Nc{p), is given by 

/+0O 
dzn(^/^T^) . (CI) 
-oo 

The mass within radius r, M(r), is given by 

M^{r) = An dnn'^ p{n) . (C2) 
Jo 

A useful quantity is the average density p. This can be defined in a number of ways, but is 
chosen here to be in terms of the half-width-at-half-maximum projected radius of the clump, 
Pi/2, defined such that A^c(Pi/2) = 0.5 Nc{p — 0). Then p is 

3 ^1/2 

The self-potential due to gravity is 

/■oo pT^ 7?'^ 

w^-igtt'^g dnn^p{n) dW —p{n') . (C4) 

Jo JQ ^ 

If Nc is defined as the column density through the clump center — i.e., Nc = Ndp = 0) — 
and Mc is defined as the total clump gas mass — i.e., Mc = Mc{r — oo), then we can write 
these and W as 

= 2kNnpi/2 , (C5) 
Me = kMPPl/2 , (C6) 

W = -kw . (C7) 

Pl/2 
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The constants k^-, kM, and kw are given by the precise functional form of p(r). Specifically, 
these constants are determined by comparing equations fICip . flC2p . and flC4p with fICSp . 
16^1 . and flC7p . respectively. 



The simple form of the Virial theorem applied in Section 12.21 and the full-width-at-half- 
maximum, di/2 = 2pi/2, are now used to derive Avc: 

Avc = k^ n5 di/2 , (C8) 

where 



kjj 



^ kw kM /umjjjG 



(C9) 



where the expressions ( IC6P and flC7p andp = fim^^ n was also used (also see Section [2^2]) . 

In addition to the above equations, we need the appropriate equation for tq. Starting 
with equation fl^ and plugging in fICSp and fICSp yields. 



which is a more general form of fl22|) . Equation fl22|) assumed a sphere of uniform density, 
for which k^ = 1, and ( IClOP reduces to ( l22l) . 

Now we can combine the expressions developed here with ( l25i) . ( l26l) . and (ITOi) to obtain, 

= (27r)^(^'^)CTA;;4U;^A;r'4~'T;^^"'n^^'"'^ • (Cll) 

This is a more general version of equation fl28|) . which assumed a uniform- density sphere. In 
fact, equation (IClip is more general than for just spherical clumps: as long as Avc, Nc{0), 
To, and A(ro) can be represented by expressions (ICSp . (ICSp . (IClOp . and (ITU]) , respectively, 
the expression for the X-factor can have the specific form given by ( IClip . 

A few examples of spherically symmetric clumps are examined in the following subsec- 
tions. 



C.l. Uniform-Density Sphere 

This case is also called the "Hard Sphere" in Section 13.11 the tables, and in MSH84. 
Because the hard sphere has a well-defined edge with a well-defined radius, R, this R is 
used in place of pi/2. The density, p(r), is constant and equal to both pc and p. Similarly 
n{r) = ric and n{r) = li for all r from zero up to R. Accordingly, it is trivial to show that 
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= kw = |, and kjy = 1. It then follows that expressions f|T6l) . fl22l) . and fl28l) result 
from ( 1C9I) . ( ICIOI) . and ( ICllI) . respectively. The central sightline optical depth, tq, is then, 

To = 1.18 X lO^n^-^T-i-^^ , (C12) 

for Tj^ = 10 to 20 K. These values along with the kA and e values given in Subsection 13.1.21 
yield, 

Xj(X2o) = 0.488 CrT-o-^6n°-^=^ , (C13) 
where this X-factor is in units of X20 or 10^^ H2 molecules ■ cm~^ ■ (K ■ km ■ s~^)~^. 
Specific numerical results of these calculations are listed in Tables [1] through [51 



C.2. Gaussian Sphere 

The optical depth profile, t{p), for the Gaussian sphere is given by ([9]). For an isother- 
mal clump in LTE, the volume density, p(r), that would give this r(p) is Gaussian: 

p{r) = poexp (^-^^ , (C14) 

where po is the density at the center of the sphere and cfr is the (l/e)-folding radius for the 
density variation. The analogous expression exists for n{r) and placing this in (IGip results 
in ^ 

Nc{p) = no V^exp (^-^^ . (C15) 
so that Nc = no ^Jii a^- From (1C15P we see that = Pi/2/Vl'n 2. The mass as a function of 



r IS 



——— = erf [ — ] — expl -] , C16 

where, 

Me(oo) = vrtpo^' • (C17) 

The erf is the Gaussian error function. Replacing r with pi/2 in ( IC16P yields MrXpi/2) I Mc{oo) 
0.2878. Placing this into (1C3I) gives usp = 0.6630 po- Accordingly, kN = 1.606, kM = 14.55. 
Using expressions (lG4p and (]G7p gives kw = \/{ln2)/2Ti = 0.3321, which in turn gives 
k^ = 3.420 X lO^^cgs. The central sightline optical depth, tq, is accordingly, 

To = 1.38 X lO^n^-^T^^-^^ , (C18) 
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for Tj^ = 10 to 20 K. These together with the fc^ and e values given in Subsection 13. 1 .31 yield. 

Xj(X2o) = 0.204 CTT-o-38nO-32 . (C19) 

Specific numerical results of these calculations are hsted in Tables [1] through [51 



C.3. Squared-Lorentzian Sphere 



This case was discussed in MSH84. The optical depth profile is given by 



rip) = ro(0) 



1 + 



p 



(C20) 



where the tq is the projected radius at which r falls to 0.25 tq. It also gives the effective 
area, where a^^^ = tttq. For an isothermal clump in LTE, it follows that 



Ncip) = iVe(O) 



1 + 



p 



2 \ H -2 



(C21) 



It is easy to show that tq = (v^ — l)~°'^Pi/2- The volume density, p(r), that corresponds to 
the surface density of equation (]C2ip is 

-2.5 



Pir) = po 



1+ -7 



(C22) 



This is easily demonstrated by substituting the expression analogous to flC22l) for nc{r) into 
IT]). Doing so gives, 

(C23) 



4 

Nc = - no ro 



Substituting flC22|) into results in, 

MJr) 



MJoo) 



+ 1 



-1.5 



where, 



Po rl 



Applying flC3l) yields p = 2~° ''^po- Consequently, 

kN = 

= 1.742 



4 / 



(C24) 



(C25) 



(C26) 



and 
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3 

26.43 



-1.5 



Plugging (IC22p into (IC4p and doing some work yields 

37r / ^ \0.5 



k 



w 



32 
0.1896 



(72-1 



From these the following is obtained: 

kv — 



IT 



(20-^^) /im,,G 



96 (72- 1) 
= 3.481 X IQ-^^cgs . 

The above results give the central sightline optical depth: 

2 -0.5 rp-1.75 



0.5 



To = 1.46 X 10^n"-°T 



K 



(C27) 



(C28) 



(C29) 



(C30) 



for Tj^ = 10 to 20 K. All of the above results along with the kA and e values given in 
Subsection 13.1.31 give 

X/(X2o) = 0.0756 CTT-o-oo3nO-22 . (C31) 
Specific numerical results of these calculations are listed in Tables [1] through [51 



C.4. p cx r ^ Sphere 

The adopted sphere density, p(r), is flat for r < tq and falls like outside of this and 
out to r = ri: 



p(r) = po for r < ro 

= Po (^) for ro<r <ri 

It is easy to show that the mass is given by 



An 



M,{n) = —porU3^ - 2 
J \ '"o 



(C32) 

(C33) 



(C34) 
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Notice that the mass, Mc(ri), would diverge if ri oo, so a maximum radius, ri must be 
specified. The usual quantities (e.g., for and kj^) must now be determined in terms of 
Mc(ri) instead of Mc{oo). Employing expressions (ICip and the number density analogs of 
l \Cm and flU33D yields 

2 norg 



Ncip) 



P 



P 



ri 



arccos — — arccos — 



P 



ro 



2-1 I P 

IT noVr^p arccos — 



for p < ro (C35) 
for rQ < p < ri, (C36) 



where N^{0) = 2noro (2 - ^). 



Determining the self-potential, W, and applying the simple form of the Virial theorem 
leads to 



'-^ - I In Mi 

Avl = AnGpoTl ^ 



14 
15 



- 2 

ro 



For simplicity, ^ 1 is assumed and 



An 

— G/im^ 



0.5 



„0.5 „ 



where po = uq fi m^^ was used. It can be shown that 



^^0 = - Pl/2 

7T 



and 



(C37) 



(C38) 



(C39) 



IT 



Po 



4(37r - 4) 
1.429 p 



(C40) 



Expressions (]C39|) and flC40p again assume that ^ ^ 1. With these expressions in hand, 
expressions for kN and ky follow: 



TT 



37r - 4 
1.819 



(C41) 



and 



kv 



3{3n - 4) 
4.20 X 10"^^ cgs 



G fim^. 



0.5 



(C42) 
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Computing Xf requires specifying ti/tq. The particular exanaple d iscussed in Sec- 
tion 14.7.11 is of the Rosette Molecular Cloud, where IWilliams et al.l (119951 ) find ^ = 5 for 
many of the clumps. Simple numerical integration yields e = 0.30 and = 1-80 for this 



value of — . Therefore, 



20 J 



0.22 Ct T, 



-0.47 -0.35 



(C43) 



Because the r ^ sphere must have an outer radius, the average density can also be 
defined as over the entire volume — i.e., n^. Doing this changes kjq and ky-. 



kN = 2- 



and 



L9 



G fim^ 



0.5 



(C44) 



(C45) 



3.18 X 10"^^ cgs 



For the clumps of the Rosette Molecular Cloud, again using the numbers of the previous 
paragraph, we have, 

Xf{X2o) = 0.88 Ct T-°-^^ nf ^ (C46) 



D. Treatment of Cylindrically Symmetric Clumps 

D.l. Viewed Perpendicularly to the Axis of Symmetry: The Side-on Case 

Here the mass density p is adopted to be a function of the radius from an axis of sym- 
metry, r, within each cylindrical clump: p = p{r). This is also true for the number density, 
n(r). Again, the x-y plane is chosen to be perpendicular to the observer's sightline and z 
is measured along the sightline. The x-axis is selected to be along the cylinder's symmetry 
axis and is also perpendicular to the observer's sightline. The y-axis is perpendicular to 
both the observer's sightline and the symmetry axis; the y is the projected distance from the 
cylinder's central axis. The radius, r, is then y"^ + z'^. The cylinder's length is h. The clump 
column density as a function of projected distance from the central axis, Nc{y), is given by 

/ + 00 
dz n{^y^ + z^) . (Dl) 
-oo 

The mass within radius r, M(r), is given by 

Mc(r) = 2n h [ dn 7^p(7^) . (D2) 
Jo 
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The half-width-at-half-maximum distance from the symmetry axis, yi/2i is defined by Nc{yi/2) 
0.5 Nc{y = 0). The average density,/), is defined analogously to that of expression flC3p : 

Mc{r = 2/1/2 

The self-potential due to gravity is 



W = -kF—^ , (D4) 

where fci? = 1 in the limit h/di/2 ^ 1 (and di/2 = 2yi/2). Equation (]D4p in this limit is 
found from solving the Poisson equation in cylindrical coordinates. The parameter /cj? is a 
correction factor for when h/di/2 ^ 1 is false; a crude numerical analysis suggests that for 
h/di/2 > ^, kp < 1.4. A similar value is found from a very crude analytical approach. If 
the cylinder is of uniform density with radius, R, and has length, h = 2R, then the cylinder 
very roughly approximates a uniform density sphere of radius, R. From Appendix IC.ll we 
know that kw = 0.6. Placing h = 2R into flD4l) and comparing with the uniform sphere, we 
find that kp ~ 1.2. In any event, the final dependence of the relevant quantities on kp will 
be weak and setting this to unity will be sufficient (see below). 

Many parameters have similar or identical expressions to those in Appendix O The 
parameter /cat is defined similarly to that in equation (1C5I) . but with yi/2 in place of pi/2- 
The kM is defined in terms of a cylindrical variation of equation ( 1C6I) : 

M, = kMphyl/^ . (D5) 

The parameter kw does not apply in this case because W is not dependent on r, except for 
the weak implicit r-dependence of kp. The expression for Av^, is the same as fICSp and the 
expression for k^ is almost the same as flC9|) : 



k 



'V 



— kp kM A* G 



(D6) 



The equations for tq and Xj are still flC10|) and flCll|) . Given that Xf oc k^^^ and that 
ky oc k^^, Xf oc A;^^*'^ ^\ Even for a value of e as low as that in a hard sphere, the 
dependence of Xf on kp is weak: Xf oc kp'^^. Since kp is between 1.0 and 1.4, the effect on 
Xf is less than 16%. 



D.1.1. Gaussian Cylinder 
The volume density, p(r), is given by, 

p{r) = p^expl-^^ . (D7) 
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This is apparently identical to flC14l) . except that the r and ar here are defined as distances 
from the central axis rather than a central point. From flD7p . it is easy to show that, 



where, 



and. 



where. 



iVc(O) = Vvr no (T^ , 



M,(r) = MJoo) 



Mc(oo) = TT po h al 



exp 



at 



(D8) 
(D9) 
(DIO) 
(Dll) 



From the above expressions, we see that yi/2 = arVln2, p = po/(2 ln2), = ln2, 
kM = 271, yielding. 



-kpp. mjj2 G 



0.5 



(D12) 



and adopting kp = 1, 



= 3.899 X IQ-^'^cgs . 
The above results give the central sightline optical depth: 



2 -0.5 rp-1.75 



To = 1.10 X lO^n^-^T 



K 



(D13) 



for Tj^ = 10 to 20 K. These results along with the /c^ and e values given in Subsection 13.1.41 
give 

X/(X2o) = 0.309 CrT-o-^6j^°-37 _ ^^4^ 



The numerical results of the above calculations are listed in Tables [T] through O 



D.2. Viewed Along the Axis of Symmetry: The End-on Case 

The mass density is again of the form p = p{r), but with the sightline along the axis of 
symmetry. The x-y plane is still chosen to be perpendicular to the observer's sightline and 
z is measured along the sightline and along the cylinder's symmetry axis. The radius, r, is 
again the distance from the central axis, but this is now + y'^- The cylinder's length again 
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is h. The clump column density as a function of projected distance from the central axis, 
Nc{r), is given by 

Nc{r) = n,(r) h , (D15) 
which is because nc{r) has no z-dependence. The HWHM radius, ri/2, is then found from 



n, 



.(ri/2) = 0.5 n,(0) . (D16) 



The expressions for Mc(r), kM,p, and are analogous to those for the side-on case — i.e., 
(ID2I) . ( ]D5I) . ( 1D3I) . and ( 1C5I) — but with ri/2 in place of yi/2 or pi/2. Combining this last 
with (IDTSD yields, 

fciv = 3 — , (D17) 

where (ii/2 = 2ri/2- The above states that oc h/di/2; i.e., /cat is proportional to the 
cylinder aspect ratio. 

Other expressions are identical to those mentioned previously: k^ is still given by flD6p . 
Avc by (ICSp . To by (IClOp . and by (IClip . The dependence of Xf on kp is the same as 
in the side-on case. The dependence of Xf on k^ implies there is now a dependence on the 
cylinder aspect ratio: Xj oc {h/di/2Y^~''^ ■ 



D.2.1. Gaussian Cylinder 

The Gaussian cylinder has been treated in the side-on case, so there are similarities in 
this end-on case. The relationship between ri/2 and ar is the same as that between 1/1/2 and 
ar- The relationship between p and po is unchanged, ku and k^ are also unchanged. One 
important change is the expression for fc^r: 

h 

kN = 2{ln2)—- . (D18) 
"1/2 

These results give tq: 

To = 1.04 X lO^nO-^T-i-^^ f-^^) , (D19) 

for Tj^ = 10 to 20 K. Another difference is that the end-on Gaussian cylinder looks like the 
Gaussian sphere. Using the kA and e for the Gaussian sphere results in, 

XfiX2o) = 0.171 CrT-'-'^n'-'' l^-^J . (D20) 
The numerical results of the above calculations are listed in Tables [T] through [5l 
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D.3. Collapsing, Magnetized Filament 



Tilley &: Pudrita (120031 ) examined the case of a constant toroidal flux-to-mass ratio 



in a collapsing cylindrical cloud. Here we estimate Xf for the side-on case only. Their 
equation (35) can be rearranged to give the velocity width: 



Av, = 2.89 X 10"^^ A 



■^frag kmax CQS 



(D21) 



where \frag is the frag mentation wavelength and kmax is a dimensionless wavenumber (see 
Tilley fc Pudrita l2003l ). They find that the density goes like for a strong magnetic field 



and r ^ for a weak magnetic field. So, assuming a roughly constant density inside radius, 

p{r) = po for r <ro (D22) 



Po 
Po 



for r > To, 



(D23) 



where a = 2 or 4. Combining equation (]D1|) with the number density analogs of (1D22P and 

mm, 



Nciy) 



oo 

j dz [z^ + y^) "'^^ for y <ro 



2 riQ^r^ + 2 nor" 

oo 


Combining equation flD2p with flD24p and flD25p results in 

Mcir) =TThrlpo + 2TThpor^ / dnn^^~"^ 



(D24) 



for y > ro. (D25) 



(D26) 



Now we will examine the a = 2 and a = 4 cases separately. 



D.3.1. a 



For a = 2 the mass diverges when r — > oo, so we must specify an outer radius, ri. The 
upper limits in equations (]D24p and (1D25P are ri instead of oo. But, for simplicity, we will 
assume that ri ^ tq. Equations (1D24I) and (1D25I) become 

2 norg 



Nc{y) 



y 



arccos | — 
ri 



arccos | — 
So 



2-1 / 2/ 

TT n^rr^y arccos 



fory<ro (D27) 
/or ro < y < ri.(D28) 
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It is easy to show that Nc{0) = AuoVq and, accordingly, that yi/2 = f tq, where it is assumed 
that ri ^ yi/2. Equation ( 1D26I) becomes 



Mc(r) = 7ihrlpo 



1 + 2ln 



(D29) 



The maximum radius, ri, must be specified and the total mass, Mc(ri), is dependent on the 
ratio of the outer-to-inner radii of the region, ti/tq. Expression (]D3I) now gives usp: 

2 



P 



I + 2ln 



7r 



Po 



(D30) 



0.771po 



With this information in hand (and assuming tq/ti to be small) it is easy to demonstrate 
that 



k 



71 



N 



1 + 2/n(f) 
1.651 



(D31) 



We can also determine ^m, but it is only useful for estimating Avc for virialized gas. Also, 
it is normally defined in terms of Mc(oo). Nevertheless, for completeness, it is given here for 
M,(ri): 

1 + 2/n(r^) 



k 



M 



IT 



1 + 2ln (f ) 



Expression (lD2ip is rewritten as 

AVe 

where, 



U f,0.5 \ J, 

fhy li ^frag i^max i 



(D32) 



(D33) 



k^ = 3.294 X 10"^^ cgs 
Starting with equation (!2TI) and plugging in ( ]D33l) and Nc{0) expressed in terms of k^ yields. 



To 



N 



^0.5 

K 



dl/2 



frag 



(D34) 



where, again, di/2 = 2?/i/2. 

Now we can combine the expressions developed here with fl25l) . f l26|l . and ffTOl) to obtain. 



Xf = i2n)^2^'-'^CTk^'krkl-'k]i'Tl^''n^^'- 



d 



1/2 



^frag kmax 



l-e 



(D35) 
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Except for the extra factor of [(ii/2/(A/rag kmax)Y~'^, this is identical to the more-or-less 
general expression for Xf, i.e., equation (IClip . Evaluating flD35p numerically depends on 
the values of and e, which depend on ri/rg. Here we try two values for this ratio: 
ri/ro = 10 and ri/ro = 1000. For ti/tq = 10, /c^ = 1.74 and e = 0.27 and, 

Xf{X2o) = 0.97 CtT:"^-^-"-^^" ' ''^ 



^-0.53 ^0.37 

K 



/2 



A 



frag 



d 



1/2 



frag 



0.31 



For ri/ro = 1000, kA = 1.46, e = 0.69 and, 

Xf{X2o) = 0.0729 CTT+o-2inOi6 

in which kmax = 0.2 was adopted for both cases (see lTilley fc Pudritdl2003l ). 



(D36) 



(D37) 



D.3.2. a = 4: 



Equations (1D24I) and ( 1D25P become 



Ar,(|/) = 2 rio^Jrl-y^ + ... 



+ no To ?/ ^ 



^ 4-3 

2 "'O'^ol/ 



arcsm 



1 - 



for y<ro (D38) 
fory>ro. (D39) 



It is easy to show that A^c(O) = | ro and, accordingly, that yi/ 
becomes 

^2 



'371^ 



MJoo) 1 - 



2rl 



2 = l"8";^ '^o- Equation ( ]D26[) 

(D40) 



where, 

Mc(oo) = 2 7r/ir^po . (D41) 

Unlike the « = 2 case, the a = 4 case has finite mass, even for r/ro cxo. Expression (1D3P 
gives us p\ 



P = '3^ 



= 0.495po 
It is now easy to demonstrate that 



1 - 1 {^^ 

2 V37r 



Po 



(D42) 



2.552 



For completeness, kM is also given: 



1 - if^)^ 

2 \3tvJ 



(D44) 



With kN known, we have Avc given by ( 1D33I) and 



4.114 X 10 



16 



cgs 



(D45) 



•V 



The equation for Xf is (ID35p . Numerically integrating ([2]) and ([3]) gives fc^ = 1-15 and 
e = 0.37. Adopting kmax = 0.2 and substituting the values for fc^ and e into ( 1D35I) results in 
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Fig. 1. — The panels above show the velocity widths, Af (^^CO), of the ^^CO J = 1 — > 
line for the Orion A and B molecular clouds versus the molecular hydrogen column den- 
sities, N(H2). The velocity widths are effective velocity widths given by the ratio of the 
velocity-integrated radiation temperature, I(CO), divided by the peak radiation temper- 
ature of the line. The N(H2) valu es for the u pper panel are those determined from the 
one-component, non-LTE models of IWalll (120061 ). The N(H2) values for the lower panel are 
the two-component, two-subsample, non-LTE models of IWalll (120061 ). The sample of points 
are those for which the intensities are 5-cr or more for the 140 fim and 240 /im continuum, the 
^^CO J = 1 — *^ line, and ^^CO J = 1 ^ line. All these maps were convolved to 1-degree 
resolution 



see 



Walll 120061 . for details). 




Fig. 2. — The effective optical depth of a clump, A(ro), after averaging over its projected 
area is plotted against the optical depth through the clump center, tq. The solid curve shows 
A(to) versus Tq for a cylindrical clump viewed perpendicularly to the symmetry axis. The 
optical depth profile across the cylinder, from the central axis towards the edges, is Gaussian. 
The dashed curve shows the corresponding curve for a spherical clump. The optical depth 
profile from the sphere center towards the edges is also Gaussian. The Gaussian spherical 
clump case was also treated and plotted in MSH84. 
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Fig. 3. — Plots of the CO J = 1 — > line radiation temperature, Tj^, normalized to its 
source function, J'^i' T^), versus the ^^CO J = 1 — derived H2 column density. Both plots 
are reproduced from [Wall (120061 ). The curves are of the form y = 1 — exp{—ax — b), w here, 
i n the ideal case, 6 = 0. The upper plot is for the LVG, one-component models of IWall 
(I2OO6I ) and the lower plot is for the LVG, two-component, two-subsample models of that 
paper. (The sour c e fun ction for the two-component models is the effective source function 



as defined in I Wall! (120061 ). ) Given that x is in units of 10 H2 cm 



(9.5±0.4)xl0-3 and 



b = (2.4 ±5.5) X 10 ^ for the upper plot and a 
for the lower plot. 



(6.4±0.2) X 10^3 and b = (7.2±0.8) x 10 



-2 
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X-Foctor (10^° cm"^ / K km s"') 



O O \J ^ O) CO 




Fig. 4. — Curves that crudely represent the X-factor's dependence on the fluffiness parameter, 
e, are depicted. The kinetic temperature, Tj^, and the average densities, represented simply 
by n, corresponding to the curves are shown. Note that the values of the X-factor for e = 
and 1 are too low by about 30%. See Section 13.2.11 for details. It must be emphasized that 
these curves are dependent on how the "average" density is defined (see Section l4.7.ip . 
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Fig. 5.— The ^^CO J = 1 ^ O/^^cq J = 1 ^ line ratio is plotted against the r^f of 
^^CO J = 1 ^ 0. The points represent positions in the Orion A and B molecular clouds 
where the peak T^C^ ^CO J = 1 ^ 0) > 3 a, and 1^(140 /im), 1^(240 /xm), and 1{^^C0) are all 
> 5 0" (see IWallll2006l ). The solid curves represent model fits that are described in Section l373l 
The reduced chi-square of these two-curve fits are given in the upper-right corner of each 
panel. The error bars are 0.5 a for a better view of the distribution of points and of the 
model curves. The upper panel has the Tg/ values derived from o ne-compon ent models and 
the lower panel has the Tg/ derived from two-component models (IWallll2006l ). 
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50 




Fig. 6.— The frequencies of values of the ^^CO J = 1 ^ O/^^cO J = 1 ^ hne ratio in the 
Orion clouds is plotted here as a histogram (sohd hne). This histogram only includes data 
where ^^CO and ^^CO observed line strengths are at or above 5-a. 
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Table 1. X- Factor Values for Different Model Clump Types 



Type of Clump'' 


kA 


e 




Xf (IG 




-7K • km 








T 


= lOK 






Tk 


= 20K 




if 2 X 10^ 


2 X lO'^ 


2 X 10 '^ 


f '= 2 X 


102 


2 X 10-^ 


2 X 10^ 


T> r 


1 





3.5 


11 


35 


1.4 




4.4 


14 


r < 1 


1 


1 


0.12 


0.12 


0.12 


0.16 




0.16 


0.16 


Hard S. 


1.7 


0.14 


1.3 


3.4 


9.1 


0.59 




1.6 


4.3 


Gaussian S. 


1.6 


0.36 


0.70 


1.5 


3.1 


0.44 




0.91 


1.9 


Sq. Lor. S. 


1.5 


0.57 


0.36 


0.58 


0.95 


0.28 




0.47 


0.76 


G. F. (side-on) 


1.5 


0.25 


0.93 


2.2 


5.2 


0.51 




1.2 


2.8 


G. F. (end-on)'i 


1.6 


0.36 


1.3 


2.8 


5.8 


0.82 




1.7 


3.6 



^These are the types of model clumps as described in Section 3: Completely Optically Thick, 
Optically Thin, Hard Sphere, Gaussian Sphere, Squared Lorentzian Sphere, Gaussian Filament 
(side-on and end-on). 

''Clump average density as defined in Section 3 and in the appendices in units of H2 molecules ■ 

'^For simplicity, a spherical geometry is adopted for computing Xf. 
Assuming a length-to-diameter ratio of 3.4. 
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Table 2. FuU-Width-at-Half-Maximum Velocity Widths of Individual Model Clumps 



Type of Clump^ 


Avc (FWHM)'' ( km 






n'== 2 X 10^ 


2 X 10^ 


2 X 10^ 


Hard S.'i 


0.46 


1.4 


4.6 


Gaussian S.'^ 


0.63 


2.0 


6.3 


Sq. Lor. S.'^ 


0.64 


2.0 


6.4 


G. Filament'' 


0.72 


2.3 


7.2 



^These are some types of model clumps as de- 
scribed in Section 3: Hard Sphere, Gaussian Sphere, 
Squared Lorentzian Sphere, Gaussian Filament. 

*^The earlier computed values have been multiplied 
by VS ln2 to convert from rms to FWHM. 

'^Clump average density as defined in Section 3 and 
in the appendices in units of H2 molecules ■ cm~^. 

•^FWHIVl diameter of 1.8 pc was adopted. 

^Length of 6.2 pc and FWHM diameter of 1.8 pc 
were adopted. 
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Table 3. Central Sightline Optical Depths, Tq, of Individual Model Clumps 



Type of Clump'^ tq 





T 


= lOK 




T 


= 20K 




11'^= 2 X 10^ 


2 X 10^ 


2 X 10^ 


11'^= 2 X 10^ 


2 X 10^ 


2 X 10^ 


Hard S. 


30 


94 


300 


8.8 


28 


88 


Gaussian S. 


34 


109 


344 


10 


32 


100 


Sq. Lor. S. 


37 


120 


370 


11 


35 


110 


G. F. (side-on) 


28 


88 


280 


8.2 


26 


82 


G. F. (end-on)'^ 


90 


280 


900 


27 


84 


270 



^These are some types of model clumps as described in Section 3: Hard Sphere, 
Gaussian Sphere, Squared Lorentzian Sphere, Gaussian Filament (side-on and end- 
on). 

^Clump average density as defined in Section 3 and in the appendices in units of 
H2 molecules • cm~^. 

^Assuming a length-to-diameter ratio of 3.4. 
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Table 4. Effective Optical Depths, A(to), of Individual Model Clumps 



Type of Clump^ kA e M'^o) 



= 10 K Tj, = 20 K 



11^= 2 X 10^ 2 X 10^ 2 X 10^ ii''^ 2 x 10^ 2 x 10^ 2 x 10^ 



Hard S. 


1.7 


0.14 


2.7 


3.2 


3.8 


2.3 


2.7 


3.2 


Gaussian S. 


1.6 


0.36 


5.7 


8.6 


13 


3.7 


5.6 


8.4 


Sq. Lor. S. 


1.5 


0.57 


12 


23 


44 


6.0 


12 


22 


G. F. (side-on) 


1.5 


0.25 


3.4 


4.6 


6.1 


2.5 


3.4 


4.5 


G. F. (end-on)^ 


1.6 


0.36 


8.0 


12 


18 


5.2 


7.9 


12 



^These are some types of model clumps as described in Section 3: Hard Sphere, Gaussian Sphere, 
Squared Lorentzian Sphere, Gaussian Filament (sidc-on and end-on). 

^Clump average density as defined in Section 3 and in the appendices in units of H2 molecules ■ 
cm~^. 

■^Assuming a length-to-diameter ratio of 3.4. 
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Table 5. Masses, Mc, of Individual Model Clumps 



lype of Clump^ Mc (Mq) 



n^= 2 X 10^ 2 X 10^ 2 X 10^ 



Hard S.^ 


3.9 


X 


10^ 


3.9 


X 


102 


3.9 


X 


10^ 


Gaussian S.'^ 


1.4 


X 


102 


1.4 


X 


10^ 


1.4 


X 


10^ 


Sq. Lor. S.^ 


2.5 


X 


102 


2.5 


X 


103 


2.5 


X 


10^ 


G. Filament*^ 


4.1 


X 


102 


4.1 


X 


10^ 


4.1 


X 


10^ 



'^These are some types of model clumps as described 
in Section 3: Hard Sphere, Gaussian Sphere, Squared 
Lorentzian Sphere, Gaussian Filament. 

^Clump average density as defined in Section 3 and 
in the appendices in units of if 2 molecules ■ cm"^. 

■^FWHM diameter of 1.8 pc was adopted. 

'^Length of 6.2 pc and FWHM diameter of 1.8 pc were 
adopted. 



